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ABSTRACT 

We  present  efficient  parallel  algorithms  for  several  basic  problems  in  computa- 
tional geometry:  convex  hulls,  Voronoi  diagrams,  detecting  line  segment  intersections, 
triangulating  simple  polygons,  minimizing  a  circumscribing  triangle,  and  recursive  data- 
structures  for  3-dimensional  queries. 
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1.  Introduction. 

Computational  geometry  addresses  algorithmic  problems  in  diverse  areas  such  as  VLSI  design, 
robotics,  and  computer  graphics.  Since  1975  there  has  been  a  wide  development  of  sequential  algorithms 
for  geometric  problems,  but  until  1985  there  was  little  published  about  developing  parallel  algonthms  for 
such  problems.  A  notable  exception  was  the  Ph.D.  research  of  Anita  Chow  (1980)  which  seems  to  have 
been  the  pioneering  work  in  the  field  but  unfortunately  only  a  portion  of  it  has  appeared  in  the  open  litera- 
ture. 

This  paper  contributes  some  parallel  algorithms  for  solving  geometric  problems: 

(1)  Convex  hulls  in  two  and  three  dimensions 

(2)  Voronoi  diagrams  and  proximity  problems 

(3)  Detecting  segment  intersections  and  triangulating  a  polygon 

(4)  Polygon  optimization  problems 

(5)  Creating  data  structures  in  two  and  three  dimensions  to  answer  some  standard  queries. 

It  is  seldom  obvious  how  to  generate  parallel  algorithms  in  this  area  since  popular  techniques  such  as 
contour  tracing,  plane  sweeping,  or  gift  wrapping  involve  an  explicitly  sequential  (iterative)  approach  (see 
[Preparaia  and  Shamos  (1985)]  for  more  details).  In  this  paper  we  exploit  data-structures  that  are  simple 
enough  to  compute  efficiently  in  parallel  while  being  powerful  enough  to  answer  queries  efficiendy.  How- 
ever, the  queries  may  be  answered  slightly  less  efficiently  than  would  be  possible  if  the  data-structures 
were  constructed  sequentially.  'Efficient'  here  means  polylogarithmic  in  parallel  time. 

All  the  algoriUims  presented  here  will  be  A^C-algorithms,  i.e.,  algorithms  executable  in  polylog  depth 
on  polynomial-size  circuits.  We  caution  that  the  circuits  here  are  slightly  different  from  those  used  in 
machine-based  complexity  theory  in  that  each  node  of  our  circuit  can  compute  an  infinite  precision  arith- 
metic operation.  Furthermore  these  algorithms  shall  all  be  described  as  if  they  were  implemented  on  a 
CREW  PRAM  -  a  concurrent-read,  exclusive-wriie  parallel  random  access  machine.  Such  a  machine  is 
conceived  as  having  a  large  number  of  processors  with  common  access  to  a  common  memory.  Any 
number  of  processors  can  read  the  same  memory  cell  simultaneously  in  0  (1)  steps,  and  any  processor  can 
write  to  a  memory  cell  in  0  (1)  steps,  but  if  two  processors  attempt  to  write  to  the  same  cell  simultaneously 
then  the  machine  enters  an  undefined  state. 

We  should  mention  here  that  Anita  Chow's  dissertation  pays  close  attention  to  the  model  of  compu- 
tation involved,  and  in  several  cases  she  provides  algorithms  for  implementation  on  two  models  of  parallel 
computation:  the  CREW  PRAM  model,  and  the  cube-connected  cycles  (CCC)  network.  As  mentioned 
before,  in  the  CREW  model,  all  the  processors  share  the  same  memory;  however,  in  the  CCC  network, 
each  processor  has  its  own  memory,  and  is  connected  to  at  most  four  other  processors.  For  details  on  the 
CCC  model,  see  [Chow  (1980)j. 

It  follows  from  [Kozen  and  Yap  (1985)]  that  most  common  computational  geometry  problems 
including  all  the  problems  considered  in  this  paper  have,  in  principle,  A^C-algorithms.  Recall  that  they 
combined  the  techniques  of  [Ben-Or,  Kozen  and  Reif  (1984)]  and  [Collins  (1975)]  to  give  parallel  algo- 
rithms for  computing  cell  decompositions  with  adjacency  information;  for  fixed  dimensions,  these  algo- 
rithms are  in  NC.   The  consequences  of  the  cell  decomposition  algorithm  for  parallel  computational 


geometry  are  explored  more  carefully  in  [Yap  (1987)].  Since  any  appeal  to  [Kozen  and  Yap  (1985)] 
involves  a  reduction  to  Tarski's  language  for  real  closed  fields,  and  the  methods  are  too  general  to  be  of 
much  practical  use,  this  furnishes  an  existence  proof  for  NC  algorithms  rather  than  furnishing  practical  NC 
algorithms. 

Many  of  the  problems  considered  here  are  known  to  have  Q(«  log  (n))  lower  bounds  in  the  algebraic 
computation  tree  model  [Ben-Or  (1983)].  Therefore  the  goal  in  our  research  in  such  cases  is  to  aim  for 
Q(log(n))  parallel  time  when  0{n)  processors  are  available.  Only  in  some  instances  do  we  present  such 
optimal  algorithms,  and  the  algorithms  may  require  a  preprocessing  step  which  involves  sorting.  While 
0{\og^{n))  sorting  networks  have  been  known  for  many  years,  to  ensure  optimality  we  may  need  to 
invoke  an  optimal  f)arallel  sorting  method  ( log  (n)  parallel  time).  The  AXS  sorting  network  achieves  these 
bounds  but  with  utterly  unacceptable  overhead  (at  the  present  state  of  knowledge:  see  [Ajtai,  Koml6s  and 
Szemer6di  (1983);  Leighton  (1984);  Cole  and  6'Dunlaing  (1986)]).  More  recendy  an  optimal  PRAM 
algorithm  has  been  described  which  has  much  lower  overhead  [Cole  (1986)].  In  our  algorithms  to  compute 
the  VcHDnoi  diagram  for  a  planar  point  set  and  convex  hull  in  three  dimensions,  we  are  able  to  avoid  sort- 
ing even  though  we  only  match  the  complexity  bounds  of  Anita  Chow. 

Also  we  shall  pay  some  attention  to  the  processor  allocation  problem:  a  PRAM  algorithm  cannot  be 
considered  completely  described  until  one  has  a  clear  description  of  how  the  various  tasks  are  to  be  allo- 
cated among  the  available  pool  of  processors.  Such  considerations  usually  become  unimportant  when  the 
time  taken  to  complete  a  parallel  step  is  logarithmic  (or  worse)  in  the  number  of  processors  involved,  since 
processor  allocation  is  usually  easy  to  accomplish  in  logarithmic  parallel  time. 

Two  standard  methods  of  parallel  algorithm  design  will  be  exploited  freely  throughout  this  paper 
computing  a  running  sum  ('parallel  prefix'),  and  list-ranking.  To  compute  a  running  sum  means  to  com- 
pute the  n  initial  sums  of  an  array  of  n  numbers.  This  is  easily  computed  in  logarithmic  parallel  time  with 
n  processors,  by  imagining  a  balanced  binary  tree  whose  leaves  are  the  array  entries:  at  each  internal  node 
the  sum  of  the  entries  of  all  its  leaf  descendants  can  be  computed  (in  an  upward  parallel  sweep),  and  from 
these  internal  sums  all  initial  sums  can  be  computed  in  a  second  phase.  List-ranking  involves  traversing  a 
linked  list  and  determining  the  rank  (distaiKe  from  the  end)  of  all  its  entries.  This  can  be  accomplished  in 
logarithmic  parallel  time  by  a  px)inter-jumping  technique;  it  is  discussed  briefly  at  the  end  of  Section  3. 

2.  Definitions  and  terms. 

As  staled  in  the  Introduction,  the  algorithms  treated  here  will  all  belong  to  the  class  NC,  but  we  shall 
express  them  as  PRAM  algorithms  running  on  a  CREW  machine.  The  notation 

NC^(/(n)) 
will  indicate  the  class  of  algorithms  running  on  a  PRAM  using /(n)  processors  and  halting  in  0(log*(/i)) 

steps.  We  are  nainly  interested  in  the  class  NC\{n),  i.e.  those  using  a  linear  number  of  processors.  Note: 
We  use  'NC*'  instead  of  'NC  to  distinguish  our  model  from  the  circuit  model  of  machine-based  complex- 
ity. See  [Yap  (1987)]. 

We  let  E^  (respectively,  E^)  denote  the  Euclidean  plane  (respectively,  Euclidean  3-space).  In  this 
paper,  we  restrict  our  attentjon  to  £^  or  £'.  Given  a  (finite  and  nonempty)  set  S  of  points  in  E^  or  £^,  the 
convex  hull  H{S)  is  the  smallest  convex  set  containing  all  points  in  5.  Sometimes  by  abuse  of  notation  we 


shall  confuse  the  convex  hull  with  the  polygon  (respectively,  polyhedron)  which  constitutes  its  boundary. 
We  use  the  notation  dH{S)  to  denote  the  polygonal  boundary  of  H  (S). 

Given  a  finite  nonempty  set  S  of  points  in  E^,  its  Voronoi  diagram  Vor(S)  is  defined  as  follows.  For 
each  point  p  in  S,  define  the  Voronoi  cell  V{p)  as  the  set  of  points  z  in  £^  such  that  z  is  as  close  to  p  as  to 
any  other  point  in  S.  Recall  that  the  Voronoi  cells  are  all  nonempty  closed  convex  sets  with  polygonal 
boundaries  and  two  cells  can  meet  only  along  their  common  boundaries  (i.e.,  the  interior  of  every  cell  is 
disjoint  from  all  other  cells).  The  Voronoi  diagram  Vor(S)  is  defined  as  the  planar  point-set  formed  from 
the  union  of  the  boundaries  of  all  the  Voronoi  cells;  it  is  a  planar  graph  with  O  (n)  edges  and  vertices,  all 
the  edges  are  straight  line-segments  (perhaps  unbounded),  and  two  cells  can  have  at  most  one  edge  in  com- 
mon. 

3.  Optimal  Parallel  Convex  Hull  in  the  Plane 

Fast  parallel  algorithms  for  the  planar  convex  hull  problem  have  been  considered  in  the  recent  litera- 
ture [Nalh  et  al.  (1981);  Chow  (1981);  Akl  (1983);  Chazelle  (1984)].  The  algorithms  typically  use  a  sim- 
ple divide-and-conquer  technique,  recursively  computing  the  convex  hull  of  n  points  by  solving  two  prob- 
lems of  half  the  size.  It  seems  that  any  such  technique  requires  Q.(\og^(n))  time.  The  main  result  in  this 
section  is  the  following: 

Theorem  1.  The  problem  of  computing  the  convex  hull  of  a  set  of  points  in  the  plane  is  in  NC^in). 

We  subsequendy  learned  that  Atallah  and  Goodrich  (1985),  and  Wagener  (1985,  1987),  have  indepen- 
dently used  a  similar  technique  to  derive  the  same  result. 

From  the  well-known  fact  that  sorting  can  be  reduced  to  the  planar  convex  hull  problem,  or  using  the 
lower  bounds  for  the  algebraic  computation  tree  model  in  [Ben-Or  (1983)],  it  follows  that  the  log(/i)  time 
is  optimal  for  n  processors.  Actually  the  algorithm  is  optimal  in  a  stronger  sense:  it  foUows  from  the  work 
of  Cook  and  Dwork  (1982)  that  0  ( log  (n))  time  is  the  best  possible  even  if  we  allow  arbitrarily  many  pro- 
cessors. 

Let  5  be  the  given  set  of  n  points.  Use  the  pair  of  points  on  the  convex  hull  //  (S)  with  maximum  and 
minimum  x-coordinate  to  partition  H{S)  into  two  parts  in  the  natural  way:  the  upper  and  the  lower  chains. 
By  symmetry,  it  is  sufficient  to  show  how  to  compute  the  upper  chain  in  0{\og{n))  time.  First  sort  S 
acccM"ding  to  the  x-coordinate  of  the  points,  using  0 (login))  time  with  n  processcrs  [Ajtai,  Komlds  and 
Szemer^  (1983),  Leighton  (1984),  Cole  (1986)].  At  this  stage,  covertical  sets  of  points  can  be  delected 
and  all  but  the  highest  from  each  set  discarded  in  0  ( log  (n))  steps.  Our  algorithm  divides  S  into  ^  sets  of 
^n  points,  recursively  computes  the  upper  convex  chain  of  each  set,  and  then  merges  all  these  chains 
together  in  0  ( log  (n))  time.  The  details  of  merging  are  non-trivial. 

Let  us  see  how  to  merge  ^n  upper  chains  in  0  ( log  (n))  steps  with  n  processors.  First  consider  two 
upper  chains  C  and  C.  Since  we  divide, S  by  x-coordinate,  this  ensures  that  the  x-coordinales  of  any  2 
chains  do  not  overlap.  Granted  that  their  edges  are  presented  in  an  array  in  sorted  order,  a  single  processor 
can  compute  in  0(log(/i))  time  the  unique  line  which  is  tangent  to  both  chains,  together  with  the  two 

points  of  tangency.  This  was  shown  by  Overmars  and  Van  Leeuwen  (1981).  Since  there  are  at  most     " 


common  tangents  to  be  computed,  n  processors  can  compute  all  such  tangents  in  time  O(log(/i)).  Note 
that  processor  allocation  raises  no  difficulty  here  since  all  the  chains  have  a  natural  sequence  (according  to 
their  x-coordinates)  and  there  is  a  natural  sequence  therefore  on  the  set  of  all  pairs  of  chains. 

Let  G  denote  the  (combinatorial)  graph  whose  vertices  are  the  points  in  S  and  whose  edges  are  the 
edges  in  the  upper  chains,  together  with  the  line-segments  defined  by  the  tangent  lines  just  discussed. 
There  are  at  most  n-1  edges  in  all  the  upper  chains,  and  they  can  be  compressed  into  contiguous  locations 


in  an  array  A  by  using  a  running-sum  technique  (in  0  ( log  (/i))  parallel  time).  The  extra 


^    <n 


- ..  edges 


may  be  appended  to  these  edges  in  the  upper  chains.  The  array  A  contains  two  directed  edges  (xj  ,0:2)  and 
{x2,xi)  for  each  undirected  edge  (j:i,ji:2).  Using  an  optimal  parallel  sorting  method,  sort  these  edges 
according  to 

(slope  of  edge)  within  (earlier  x-coordinaie). 
The  slope  of  an  edge  is  the  angle  it  makes  with  the  x-axis,  measured  from  -7:/2  to  3n/2  (with  k/2 
ccOTesponding  to  the  vertical  upwards  direction).  Thus  the  array  is  partitioned  into  adjacency  lists  for  the 
vertices  of  G. 

From  G  we  can  identify  vertices  of  the  upper  chain  of  S  using  the  following  observation:  Consider 
any  point  v  in  5,  which  is  neither  leftmost  nor  rightmost  in  S.  If  v  belongs  to  the  upper  chain  of  S,  let  u  and 
w  be  the  vertices  immediately  before  and  after  v  along  the  chain.  Clearly  the  slopes  of  the  edges  (v,u)  and 
(v,w)  are  adjacent  among  ail  the  edges  out  of  v,  and  the  angle  uvw  (measured  clockwise  from  u  to  w)  is 
>  7t.  Clearly,  for  each  vertex  v  there  can  be  at  most  one  pair  of  edges  (v,w)  and  (v.w)  with  this  property, 
and  such  vertices  v  can  be  detected  in  constant  lime  with  one  processor  per  edge  in  the  array  A.  This  cri- 
terion identifies  all  edges  and  vertices  on  the  upper  chain,  and  perhaps  some  others,  except  far  the  first  and 
last  edges  on  the  upper  chain:  these  are,  of  course,  the  edges  with  maximum  and  minimum  slopes  respec- 
tively leaving  the  leftmost  and  rightmost  vertices.  (Since  S  is  presorted  by  x-coordinate  it  is  easy  to  iden- 
tify these  two  vertices.) 

For  all  these  triples  u,v,w  which  indicate  possible  vertices  of  the  upper  chain,  let  us  call  v  a  'candi- 
date' and  w  its  'successor.'  If  v  is  not  a  candidate  its  successor  is  undefined,  and  the  leftmost  vertex  has  as 
its  successor  the  next  vertex  on  the  upper  chain.  As  indicated  above  this  is  easily  identified.  This  succes- 
sor relation  defines  a  forest  of  trees,  in  which  the  upper  chain  is  represented  by  a  braiKh  connecting  the 
leftmost  vertex  to  the  rightmost  This  branch  can  be  identified  in  logarithmic  parallel  time,  one  processor 
per  vertex,  using  the  parallel  ranking  technique  described  below. 

To  complete  our  inductive  step,  we  now  must  collect  the  vertices  of  the  upper  chain  of  S  into  con- 
secutive locations  in  an  array.  This  is  easily  accomplished  using  parallel  prefix,  i.e.,  computing  a  running 
sum  on  an  array  C(v)  whose  entry  is  1  if  v  is  in  the  upper  chain  and  0  otherwise.  This  completes  the  merge 
process  of  our  algorithm. 

Since  the  ranking  problem  is  so  basic  and  will  be  used  again,  it  bears  repeating  here.  Here  we  con- 
sider the  successor  relation  as  defining  a  forest  of  trees  directed  from  child  to  parent,  and  the  rank  of  a  node 
in  this  forest  is  its  distance  from  its  root  ancestor.  In  view  of  a  later  application,  in  addition  to  the  rank 
information,  we  compute  an  additional  0  ( log  (n))  pointers  for  each  vertex: 


(*)      In  log(n)  parallel  time,  each  node  can  leam  its  rank  and  also  for  each  i  <  log(n),  it  also  holds  a 

pointer  to  its  2'-th  successor,  if  it  has  such  a  successor. 
This  is  done  in  two  phases,  each  phase  taking  log  (n)  parallel  time.  In  the  first  phase,  by  the  well-known 
doubling  technique,  every  vertex  determines  for  all  j  <  log(;i)  whether  it  has  a  2'-th  successor,  and  if  so, 
which.   In  the  second  phase,  at  stage  /  (there  are  log(n)  stages)  all  vertices  with  rank  between  2'  and 
2'*'-l  leam  their  exact  rank.  It  is  easy  to  implement  all  this  on  a  CREW  PRAM  in  logarithmic  time. 

The  runtime  T{n)  of  this  convex  hull  algorithm,  for  n  points  with  n  processors,  satisfies  the  relation 

r(n)<c  log  (n>+7"(^) 
which  implies  T{n)  is  0(log(n)). 

4.  Voronoi  Diagram  and  Proximity  Problems 

One  of  the  fundamental  data  structures  of  computational  geometry  is  the  Voronoi  diagram  [Shamos 
(1977)].  The  first  optimal  sequential  algorithm  to  compute  the  Voronoi  diagram  was  given  by  Shamos  and 
Hoey  (1975)  and  it  takes  0(n  log  n)  lime.  In  this  section,  we  show 

Theorem  2.  There  is  an  NCjin)  algorithm  for  the  planar  Voronoi  diagram  of  a  set  of  points.  Further- 
more, the  algorithm  only  assumes  an  A'C  2(1)  sorting  algorithm. 

Shamos  (1977)  also  pointed  out  that  the  Voronoi  diagram  can  be  used  as  a  powerful  tool  to  give  efficient 
algorithms  for  a  wide  variety  of  other  geometric  proximity  problems.  From  the  Vc^onoi  diagram  we  can 
easily  obtain  the  closest  pair  of  sites  (points),  or  the  closest  neighbor  of  each  site,  or  the  Euclidean 
minimum  spanning  tree  of  the  n  sites,  or  the  largest  point-free  circle  with  center  inside  their  convex  hull, 
etcetera.  Efficient  reductions  for  these  problems  are  possible  in  our  parallel  computation  model,  so  our 
polylog  result  for  the  Voronoi  also  implies  polylog  algorithms  for  all  these  problems,  using  a  linear  number 
of  processors.  Note  that  there  is  a  well-known  reduction  of  the  Voronoi  diagram  of  a  set  of  points  in  the 
plane  to  the  convex  hull  of  a  set  of  points  in  3-space  [Brown  (1979);  Guibas  and  Stolfi  (1983)].  Also,  in 
her  thesis,  Anita  Chow  (1980)  shows  that  the  three-dimensional  convex  hull  can  be  computed  in  NC'2(.n), 
and  exploits  this  to  give  an  0(log^(n))  algorithm  to  compute  the  diagram.  However,  her  methods  fre- 
quently involve  optimal  time  parallel  sorting',  which  the  following  direct  method  avoids.  In  our  FOCS 
extended  abstract  we  gave  an  NCjin)  algorithm  which  avoids  optimal  time  parallel  sort.  In  the  present 
paper,  we  improve  this  to/^C2(n)  by  exploiting  a  simple  observation  suggested  by  Goodrich. 

The  present  algorithm  for  computing  the  Voronoi  diagram  is  modelled  after  the  sequential  method 
first  presented  by  Shamos  and  Hoey  [1975].  This  is  a  divide-and-conquer  method  in  which  two  Voronoi 
diagrams  for  the  respective  sets  P  and  Q  of  points  (where  these  sets  are  vertically  separated,  i.e.,  there  is  a 
vertical  line  L  which  has  all  of  f  to  its  left  and  all  of  Q  to  its  right)  are  merged  in  Unear  lime  to  a  diagram 
for  P  uQ.  The  merging  is  defined  by  a  polygonal  path  (the  'contour'  or  'seam')  separating  P  from  Q;  the 
process  is  akin  to  tracing  a  line  and  seems  inherently  sequential.  The  contour  separates  the  plane  into  the 
set  of  points  closer  to  P  than  to  Q  (i.e.,  the  set  of  points  x  whose  closest  point  or  points  in  5  are  all  in  /")  and 


'  The  mmime  bound  given  in  [Chow  (1980)]  diffen  from  the  bound  quoted  here  by  an  additional  factor  of  log  log  (n). 
These  discrepancies  ariM  wherever  sorting  is  involved.  The  0{log(rt)log  log(/t))  parallel  sorting  algorithm  in  [Valiant 
(1975)1  has  since  been  superseded  by  theoretically  optimal  algorithms. 


those  closer  to  Q  than  to  P. 

The  following  illustrates  the  contour  between  two  hypothetical  sets  P  and  Q: 


contour 


Figure  1 .  Merging  by  contour-tracing 
Our  goal  is  to  parallelize  this  contour  tracing.  Let  us  introduce  some  suitable  terminology.  The  con- 
vex hull  H  (S)  has  already  been  defined.  The  normals  of  S  are  defined  as  the  half-lines  extending  from  the 
comers  of  H{S)  and  perpendicular  to  the  sides  incident  to  those  comers:  two  normals  extend  fit>m  each 
cwner.  These  normals  partition  the  complement  of  //(S),  E^~H{S),  into  sectors,  which  are  alternately 
slices  (V-shaped)  and  snips  (bounded  by  a  side  oi  H(S)  and  the  two  incident  normals).  The  following  pic- 
ture clarifies  this  terminology: 


x-/ 


Figure  2.  Normals  and  sectors  (slices  and  strips) 

The  points  of  intersection  of  Vor{S)  with  the  convex  polygon  dH  (5)  we  call  the  cuipoints  (figure  3).  Let  U 
be  the  strip  defined  by  an  edge  e  of  the  convex  hull.  The  set  t/  n  Vor{S)  we  call  the  e-promoniory.  We 
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will  assume  thai  S  is  in  general  position  (no  three  points  collinear,  no  four  concyclic).  Together  with  the 
convexity  of  the  Voronoi  cells,  this  implies  that  the  promontory  has  the  natural  structure  of  a  binary  tree 
with  leaves  at  the  cutpoints  and  such  that  the  parent  of  a  node  is  always  more  distant  from  e.  Beyond  the 
root  of  this  binary  tree  extends  an  infinite  ray  that  bisects  the  strip. 


Figure  3.  Promontory:  cutpoints  are  indicated  by  (•) 
Point  location  inside  strips  is  easily  done  in  0  (log^(/:))  serial  time.  We  now  modify  a  suggestion  of 

Goodrich  to  show  that  O  ( log  (n))  serial  time  suffices. 

Lemma  3.  Suppose  a  strip  has  an  e-promontory  with  k  vertices  and  edges.  There  is  ak  processor 
0(log()t))  time  parallel  algorithm  for  constructing  a  data-structure  for  the  plane  partition  defined  by  the 
e-promontory  so  that  in  0  ( log  {k))  time  a  single  processor  can  answer  a  point  location  query. 

Proof.  The  e-promontory  is  a  complete  binary  tree  laid  out  in  the  plane.  We  may  assume  that  each  leaf  is 
on  the  j-axis,  and  for  every  node  v  all  its  descendants  have  smaller  >'-coordinate.  For  any  node  v  of  the 
promontory,  let  Ty  denote  the  subtree  of  the  c-promontory  rooted  at  v.  Each  internal  node  v  of  the  e- 
promonlory  tree  may  be  associated  with  a  triangular  region  R^  with  v  as  one  comer  of  the  triangle  and  the 
two  other  comers  the  extreme  cutpoints  below  the  subtree  at  v.  By  convexity  of  Voronoi  cells,  the  entire 
subtree  T,  is  contained  in  /?,.  Suppose  p  is  a  query  point  inside  the  strip.  Let  r  be  the  root  of  the  e- 
promontory.  The  interesting  case  is  when  p  lies  inside  R,,  and  this  can  be  checked  in  constant  time.  We 
can  now  do  a  binary  search  inside  R,  by  locating  a  vertex  v  such  that  T,  contains  between  1/3  and  2/3  of 
the  vertices  of  T,..  If  p  is  inside  7,,  we  can  recursively  do  point  location  for  p  inside  this  subtree.  Other- 
wise, let  we  next  search  inside  a  modified  version  of  T/.  the  modified  tree  T^  consists  of  replacing  T,  by 
three  vertices  corresponding  to  the  comers  of  the  triangle  /?,,  together  with  the  two  edges  from  v  to  the 
other  two  comers.  This  recursive  search  can  be  carried  out  in  0  ( log  {k))  serial  time. 

It  remains  to  construct  a  search  structure  D,  for  the  tree  T,.  The  search  stmcture  D^  consists  of  a 
binary  tree  whose  internal  nodes  are  associated  with  triangular  regions  of  the  form  /?,  which  we  use  for 
comparisons:  if  we  are  at  some  internal  node  z  of  D„  we  next  go  to  the  right  (resp.  left)  child  of  z  if  query 
point  p  lies  inside  (resp.  outside)  the  triangular  region  of  z.  To  construct  D^  associated  with  T^,  we  do  an 
important  preprocessing  step  where  for  each  node  v  we  determine  its  preorder  rank  rank(v)  in  the  tree  and 
from  this,  we  easily  deduce  maxrank(v),  the  maximum  rank  of  its  descendants.  This  can  be  done  using  k 
processors  and  0(log(Jk))  time  [Taijan  and  Vishkin  (1985)].^  We  will  deduce  two  important  pieces  of 


^  Their  CRCW  parallel  RAM  algorithm  can  be  adapted  for  our  CREW  version. 


information  from  these  computed  values.  First,  the  number  desc  (v)  of  descendants  of  a  node  v  is  given  by 
\+maxrank{vy-rank{v).  Second,  we  can  decide  in  constant  time  if  any  given  vertex  u  is  a  descendant  of 
another  v.  More  precisely,  u  is  a  descendant  of  v  iff 

rank  (v)  <  rank  (u)  <  maxrank  (v). 
Inductively,  suppose  we  have  to  construct  a  search  structure  D^  for  the  subtree  Tr  rooted  at  some  node  r. 
We  assume  that  desc  (v)  is  precompuled  for  each  v  in  T,.  In  particular,  let  k  =  desc{r)  be  the  number  of 
nodes  in  T,.  Assigning  one  processor  to  each  vertex  of  T,.  in  0  (1)  steps,  we  can  locate  the  deepest  node  of 
T,  with  at  least  2jt/3  descendants,  and  choose  the  node  vq  to  be  that  one  of  its  children  with  the  larger 
number  of  descendants  (breaking  ties  arbitrarily).  It  is  easy  to  see  that  vq  will  have  between  i/3  and  2fc/3 
descendants.  The  root  of  D,  will  be  associated  with  the  region  R^.  The  right  subtree  of  (the  root  oO  Dr 
will  be  recursively  constructed  from  7,^:  note  that  the  precomputed  information  desc  (v)  icx  each  node  of 
r,g  is  available  with  no  further  work.  Let  T,.  be  the  modified  tree  obtained  by  pruning  r,g  as  described 
above.  We  want  to  recursively  construct  the  left  subtree  of  D,  from  T^.  But  first  we  must  revise  the 
precomputed  information  desc{v)  for  each  node  for  T^.  To  do  this,  we  only  need  to  update  the  value  of 
desciy)  by  subtracting  desc{vQ}-2  from  desc{v)  for  all  v  that  are  ancestors  of  vg.  But  this  is  easy  since 
with  one  processor  per  vertex,  each  vertex  can  in  constant  time  deduce  if  it  is  an  ancestor  of  vg  (use  the 
preorder  numbering  of  nodes  as  shown  above).  Note  that  to  represent  the  partition  of  the  nodes  of  T,  into 
two  sets  (belonging  to  T^^  and  to  T,  respectively),  it  is  enough  to  have  each  node  of  T,  keep  a  pointer  to 
one  of  the  children  of  the  root  r  in  D^.  The  region  R^  can  be  written  into  a  location  associated  with  r. 
Q£.D. 

The  above  data-structure  will  help  us  to  form  the  contour  efficiently.  Let  us  introduce  some  simple 
notation.  Recall  that  a  recursive  step  in  the  algorithm  will  involve  separating  a  given  set  5  by  a  vertical  line 
L  into  two  sets  /"uQ  (where  ?  is  to  the  left  of  L  and  Q\o'\\s  right),  recursively  constructing  Vor{P)  and 
Vor{Q)  separately  and  merging  the  two  along  the  contour.  We  shall  speak  of  the  edges  of  Vor{P)  and 
Vor{Q)  as  /'/'-edges  and  QQ-odgts.  In  the  combined  Voronoi  diagram  the  contour  edges  can  be  called  PQ 
edges  in  an  obvious  extension  of  this  notation,  and  the  vertices  can  be  classified  as  PPP,  PPQ,  PQQ,  and 
QQQ,  depending  on  whether  the  associated  three  points  of  5  (assuming  general  position  so  no  four  are 
cocyclic)  involve  3,  2,  1,  or  no  points  of  P. 

We  need  therefore  to  identify  the  PPQ  and  PQQ  vertices.  Instead  of  computing  them  directly,  we 
solve  the  following  simpler  problem:  to  identify  the  PP  and  QQ  edges  that  intersect  the  contour.  We  post- 
pone till  later  determining  precisely  where  they  meet  the  contour. 

To  carry  this  out,  we  need  one  more  preliminary  computation.  This  is  necessitated  by  the  fact  that  a 
Voronoi  edge  may  intersect  the  contour  twice.^  However  we  can  break  up  each  Voronoi  edge  into  at  most 
two  'semi-edges'  such  that  the  contour  meets  each  semi-edge  at  most  once:  let  e  be  any  edge  of  Vor(P) 
determined  by  the  two  points  u,  v  in  P.  If  the  horizonial  ray  emanating  leftward  fi^om  u  or  v  intersects  e  at 
some  point  x  (this  happens  with  at  most  one  of  the  rays  from  u  and  v)  then  we  split  e  into  two  semi-edges  at 
X.  If  neither  ray  intersects  e  then  the  entire  edge  e  constitutes  a  semi-edge.  Similarly,  if  e  is  an  edge  of 
^or(Q),  we  would  split  e  into  at  most  two  semi-edges  using  the  horizontal  rays  emanating  rightward  from 


'  We  «re  indcbtul  to  Hubert  Wagcner  (pnvue  communiuiion)  for  poiniing  this  out  to  us. 
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uor  V. 

Lemma  4.  Each  semi-edge  intersects  the  contour  at  most  once. 

Proof  (sketch).  Let  e  be  an  edge  separating  points  u  and  v  of  /*  in  VoriJP).  Since  any  intersection  of  the 
contour  with  e  results  in  a  PPQ  vertex,  which  will  be  an  endpoint  of  an  edge,  call  it  e',  separating  u  and  v  in 
the  combined  diagram,  we  conclude  that  the  contour  cannot  intersect  e  more  than  twice.  Suppose  that  it 
intersects  e  nvice.  Then  e'  together  with  the  contour  would  enclose  one  of  these  points,  u  say,  completely, 
and  since  the  contour  is  monotonic  in  the  ^-direction  it  is  easy  to  see  that  the  horizontal  line  extending  left- 
wards from  u  intersects  the  edge  e'.  Thus  each  intersection  point  is  on  a  different  semi-edge  of  e.  Q£.D. 

Call  any  semi-edge  that  intersects  the  contour  relevant.  Note  that  an  semi-edge  in  Vor{f)  is  relevant 
if  and  only  if  one  endpoint  is  closer  to  P  and  the  other  endpoint  is  closer  to  2.  To  decide  if  a  Voronoi  ver- 
tex (or  other  point  bounding  a  semi-edge)  v  of  Vor{P)  is  closer  to  P  than  to  Q,  we  have  two  cases,  (i)  If  v 
is  left  of  L  then  we  determine  the  slice  or  strip  oi  H{Q)  containing  v.  If  v  is  in  a  slice,  then  v  is  closest  to  a 
unique  point  q  in  Q;  we  can  decide  easily  if  v  is  closer  to  q  than  to  the  closest  points  of  P.  If  v  is  in  a  strip, 
we  reduce  the  problem  to  point  location  inside  a  strip;  once  we  have  located  v  in  the  appropriate  ceU  of 
Vor{Q)  we  can  compare  distances  as  in  the  case  of  the  slice,  (ii)  If  v  is  right  of  L  it  may  initially  appear 
that  we  need  a  general  point  location  data-structure  for  that  part  of  Vor{Q)  lying  to  the  right  of  L.  We 
show  how  to  avoid  this  using  another  idea. 

We  need  some  new  terminology.  Consider  the  set  of  points  x  to  the  right  of  L  such  that  the  circle 
centered  at  x  and  tangential  to  L  passes  through  some  point  q&Q  but  where  the  interior  of  the  circle  con- 
tains no  other  points  of  Q.  The  set  of  such  points  will  be  call  the  beachline  (of  Q).  (Alternatively,  we 
could  have  defined  the  beach  to  be  the  'contour'  separating  L  from  Q.)  Fcr  any  point  ^  not  in  L,  let  the  q- 
parabola  be  the  parabola  focused  at  q  with  directrix  L.  (See  figure  4). 


Figure  4.  The  beachline  (the  beachhead  B  is  shaded). 


Lemma  5. 


(a)  The  beachline  is  a  simple  curve,  disjoint  from  L  and  projecting  monotonically  onto  L.  Thus  the 
beachline  divides  the  half-plane  to  the  right  oft  into  two  connected  regions.  The  region  that  is  adja- 
cent to  L  is  called  the  beachhead  and  denoted  by  B.  The  beachhead  is  characterized  as  the  set  of 
those  points  in  the  right  half-plane  ofL  that  are  closer  to  L  than  to  any  point  in  Q. 

(b)  The  beachline  is  composed  of  a  finite  sequence  of  parabolic  arcs,  each  such  parabola  being  a  par- 
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tion  of  some  q-parabola,  qe  Q.  Each  such  q  is  said  to  be  near  L.  The  point  where  two  parabolic 
arcs  meet  is  called  a  transition  point. 

(c)  Let  C  be  the  contour  defined  by  the  partition  S-P  <jQ.  Then  C  lies  entirely  to  the  left  of  the  beach- 
line. 

(d)  The  Voronoi  cells  of  Vor  {Q)  that  intersect  B  are  precisely  those  belonging  to  points  q  that  are  near 
L  in  the  sense  of  (b).  Each  connected  component  in  the  restriction  of  Vor  {Q)  to  B  has  the  structure 
of  a  free  tree  embedded  in  the  plane  whose  leaves  are  on  the  boundary  ofB.at  most  one  lec^is  on  L. 
The  intersection  of  a  cell  of  Vor  {Q)  with  the  beachhead  may  have  several  connected  componerus. 
Each  such  component  is  called  a  B-cell.  The  B-cells  are  linearly  ordered  along  the  beachline. 

(e)  Let  X  be  a  transition  point.  Then  there  are  two  points  q,q'^  Q  such  that  x  lies  on  an  edge  of  Vor  (Q) 
separating  the  cells  of  q  and  q'.  Moreover,  the  two  parabolic  arcs  that  meet  at  x  are  portions  of  the 
q-  and  q' -parabolas,  respectively. 

Proof,  (a)  To  see  that  the  beachline  projects  monotonically  onto  e,  for  each  point  p  on  L,  draw  the  horizon- 
tal ray  ^  emanating  rightward  from  p.  The  beach  line  interesects  /?  at  a  unique  point  x  such  that  the  points 
on  R  closer  to  L  than  to  2  are  precisely  those  lying  between  p  and  x. 

(b)  The  points  x  of  the  beachline  where  the  circle  centered  at  x  and  tangent  to  L  is  incident  to  two  or 
more  points  of  Q  are  finitely  many  in  number.  These  define  the  transition  points.  For  non-transition  points 
X  on  the  beachline,  the  circle  at  x  touches  a  unique  q&Q.  Such  an  x  lies  on  the  parabola  focused  at  q.  Thus 
the  beachUne  is  composed  of  a  finite  number  of  parabolic  arcs  that  are  joined  together  at  transition  points. 

(c)  Let  j:  be  a  point  on  the  contour,  and  let  D  be  the  largest  disc  centered  at  x  whose  interior  does  not 
contain  any  point  of  PkjQ.  Then  D  touches  points  both  in  P  and  Q  since  x  is  on  the  contour,  therefore  its 
interior  intersects  L.  Were  x  to  the  right  of  the  beachline  it  would  follow  that  D  contained  points  of  2  in  its 
interior,  a  contradiction;  hence  j:  is  on  or  to  the  left  of  the  beachline. 

(d)  Qearly  a  point  in  2  is  near  L  if  and  only  if  its  cell  in  Vor  (Q)  intersects  B.  Since  every  point  in  Q 
is  obviously  to  the  right  of  the  beachline,  every  fl-cell  must  intersect  the  beachline  and  therefore  the  5-cells 
are  linearly  ordered  by  their  intersections  with  the  beachline.  For  the  same  reason  each  connected  com- 
ponent of  Br\Vor(Q)  is  acyclic  and  has  the  structure  of  a  (free)  tree  whose  leaves  are  on  the  boundary  of 
B.  If  such  a  component  met  the  line  L  twice  then  there  would  be  a  chain  of  edges  connecting  these  two 
points  within  B  and  therefore  there  would  exist  some  B-cells  not  meeting  the  beachline:  hence  at  most  one 
leaf  is  on  L. 

(e)  Choose  q  and  <?'  to  be  the  points  that  define  the  parabolic  arcs  that  meet  at  x.  QE.D. 

Lemma  6.  InOi  log  (n))  time  with  0  (n)  processors,  we  can  construct  the  beachUne  of  Q  and  a  suit- 
able data-structure  to  do  the  following  query:  For  any  pnint  v,  to  determine  in  0{\og{n))  serial  time, 
whether  v  lies  in  the  beachhead  B  and,  if  so,  which  B-cell  contains  v. 

Proof.  Using  one  processor  per  edge  of  Vor{Q),  we  can  compute  all  the  transition  points  (at  most  2)  on 
that  edge.  To  see  this,  let  e  be  an  edge  separating  the  Voronoi  cells  of  q  and  q'  in  Vor{Q).  Then  at  each 
point  X  (at  most  2)  where  the  (^-parabola  intersects  e,  we  get  a:  is  a  transition  point  Once  all  the  transition 
points  are  computed,  we  may  sort  them  by  their  y-coordinates  (using  list-ranking)  to  form  the  beachline  as 
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well  as  determining  all  the  B-cells. 

Next  consider  the  point  location  problem:  by  lemma  5(a),  we  can  use  binary  search  to  decide  if  a 
query  point  p  lies  in  the  beachhead  in  0  ( log  {n))  time.  It  is  slightly  more  involved  to  locate  the  fi-cell  con- 
taining p  in  the  beachhead. 

Recall  that  Vor{Q,)r-£  is  composed  of  a  collection  of  free  trees  with  leaves  on  the  beachline  and 
perhaps  one  on  L.  We  can  assign  a  root  to  each  such  free  tree:  if  it  reaches  L,  let  its  root  be  the  vertex  con- 
nected to  L;  otherwise,  let  its  root  be  the  highest  leaf  vertex,  say.  AH  the  edges  then  can  be  directed 
towards  the  root  using  a  trick  of  Tarjan  and  Vishkin's  (1985):  form  an  eulerian  directed  path  beginning  at 
the  root  and  visiting  each  edge  twice,  in  which  case  v  is  an  ancestor  of  w  if  and  only  if  v  is  first  visited 
before  w  and  last  visited  after  w,  thus  the  edges  can  be  directed  using  list  ranking.  We  can  form  regions  ^7- 
that  enclose  each  individual  tree  T:  If  T  reaches  the  line  L  then  the  region  ^7-  has  three  comers  formed  from 
the  root  of  the  tree  and  the  two  extreme  leaves  of  the  tree.  If  T  does  not  reach  L  then  let  u  and  v  be  its  top 
and  bottom  intersections  with  the  beachline.  The  line  joining  these  two  points  cannot  be  crossed  by  any 
edge  of  T  since  such  an  edge  would  have  to  cross  the  boundary  of  fi  in  L  or  outside  the  range  of  these 
points;  hence  all  of  7  is  to  the  right  of  this  line  and  we  define  Rj  as  that  region  of  B  to  the  right  of  this  line 
in  this  case.  (Unbounded  regions  raise  a  few  special  cases  which  are  easily  dealt  with.)  Each  connected 
region  of  the  beachhead  exterior  to  all  the  regions  Rj  belong  to  a  unique  fi-celL  Thus  using  binary  search 
we  can  determine  the  region  R-j-  (if  any)  that  contains  the  query  point  p.  If  p  lies  outside  all  the/?  7-  then  we 
are  done.  Otherwise,  we  use  the  same  technique  as  for  point  location  in  an  e-promontory.  Q.E.D. 

We  are  now  ready  to  give  the  overall  algorithm  in  five  steps  (A)-(F).  Step  (A)  is  done  once  and  the 
rest  are  recursively  called.  The  complexity  0  ( log  ^(n))  follows  since  steps  (B)-(F)  each  take  0  ( log  (/i)). 

(A)  First  sort  S  according  to  its  jr-coordinate  (an  0  ( log  ^(n))  sorting  method  is  enough  here).  Steps  (B)- 
(F)  involve  partitioning  S  evenly  into  F  u  2  (using  a  vertical  line  L)  and  recursively  computing 
Vor{P),  Vor{Q)  and  also  the  convex  hulls  //(?),  H(Q). 

(B)  Construct  the  beachline  for  Q  and  the  corresponding  search  structure  for  point  location  within  the 
beachhead.  Similarly,  construct  the  data  structure  for  point  location  within  each  strip  emanating 
from  the  convex  hull  of  Q.  This  takes  0  ( log  (n))  time  with  n  processors.  Do  the  same  for  P. 

(C)  For  each  point  v  bounding  a  semiedge  of  Vor  (P)  determine  whether  it  lies  closer  to  Q  than  to  P, 
using  the  two  techniques  shown  above:  if  v  lies  to  the  left  of  L,  determine  the  sector  (strip  or  slice) 
containing  v  and  if  a  strip,  do  a  point  location  to  determine  the  cell  containing  v.  If  v  lies  to  the  right 
of  L  then  we  check  if  v  lies  within  the  beachhead,  and  if  not  we  know  that  v  is  closer  to  Q  than  to  P. 
Otherwise,  we  locate  the  cell  containing  v.  In  either  case,  we  can  decide  if  v  is  closer  to  Q  than  to  P. 
This  can  be  accomplished  in  0  ( log  (n))  time  with  one  processor  per  vertex.  Repeat  the  step  with  the 
roles  of  P  and  Q  reversed. 

(D)  At  this  point  we  know  which  semiedges  of  Vor(P)  (and  by  a  symmetrical  operation,  Vor{Q))  meet 
the  contour.  Using  this  information,  we  shall  next  compute  a  linear  sequence  of  triangular  subcells 
of  cells  of  Vor{P),  whose  union  covers  the  contour  and  such  that  the  contour  visits  each  triangular 
subcell  exactly  once,  and  visits  the  cells  in  the  given  linear  order.  We  call  this  structure  the 
P -conduit.  The  contour  will  finally  be  determined  by  combining  the  P-conduit  with  the  2-conduit 
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(a  merging  process). 

It  will  simplify  things  if  all  the  geometric  constructions  are  bounded.  This  can  be  ensured  by  finding 
a  rectangle  which  contains  on  its  inside  all  points  in  S,  and  also  the  top  and  bottom  contour  vertices.  Since 
the  beachlines  on  both  sides  are  enclosed  by  the  vertical  sides  of  such  a  rectangle,  we  conclude  thai  the 
contour  can  only  cross  its  horizontal  sides,  and  no  essential  information  is  lost.  In  our  discussion  of  the 
conduits  we  shall  assume  that  all  the  cells  and  edges  are  bounded,  having  been  'clipped'  to  lie  inside  this 
rectangle.  It  is  relatively  easy  to  construct  such  a  rectangle  in  logarithmic  parallel  time.  Details  are  left  to 
the  interested  reader. 

Consider  a  relevant  point  p  of  P,  i.e.,  one  whose  cell  meets  the  contour.  The  contour  first  enters  the 
p-cell  (in  Vor{P),  perhaps  clipped)  through  an  upper  semiedge,  possibly  leaves  and  returns  by  several  pairs 
of  semiedges  (each  pair  being  two  halves  of  the  same  edge),  and  finally  leaves  the  cell  for  ever.  The  direc- 
tion of  travel  within  the  p-cell  is  clockwise  around  p  and  all  parts  of  the  contour  within  the  cell  are,  of 
course,  visible  from  p.  Let  e  and  /be  a  contiguous  pair  of  semiedges  in  this  clockwise  sequence  where  e  is 
an  'entering'  semiedge  and  /a  'leaving'  semiedge.  Consider  these  two  semiedges  together  with  the  list  of 
all  edges  bounding  the  (clipped)  p-cell  included  between  them  in  clockwise  order  around  p.  Connecting 
the  endpoints  in  this  ordered  list  of  edges  to  p  (including  e  and  f),  we  obtain  an  ordCTed  list  of  triangles  with 
the  property  that  the  contour  enters  the  first  and  leaves  the  last  (to  enter  a  different  cell)  and  enters  and 
leaves  each  intermediate  triangle  exactly  once.  This  list  of  triangles  forms  part  of  the  P-conduit  mentioned 
above.  It  is  relatively  straightforward  to  combine  all  these  lists  together  for  all  relevant  points  p,  in  the 
appropriate  linear  sequence,  by  assigning  the  appropriate  number  of  processors  to  each  point  p:  namely, 
one  to  each  edge  of  the  clipped  p-cell,  thereby  locating,  counting,  and  sequencing  the  lists  of  triangles 
around  p:  to  link  these  lists  together  in  the  correct  way  is  easily  done  using  one  processor  per  semi-edge. 
The  time  to  do  this  is  O  ( log  (/»)).  Thus  the  P-conduit  can  be  formed  as  a  linked  list  and  we  can  use  list- 
ranking  to  store  an  appropriate  representation  of  this  structure  in  an  indexed  array. 

(E)  Finally  we  construct  the  contour  by  merging  the  /'-conduit  with  the  Q-conduiL  Say  that  two  triangles 
Tp  and  Tq  from  the  respective  conduits  interact  if  the  bisector  separating  the  respective  points  p  and 
q  (apexes  of  the  respective  triangles)  meets  the  intersection  of  the  two  triangles.  Clearly,  these  trian- 
gles interact  if  and  only  if  the  contour  contains  an  edge  separating  the  two  points  and  this  edge  inter- 
sects the  common  intersection  of  both  triangles.  Each  conduit  triangle  interacts  with  at  least  one  on 
the  other  side.  Thus  the  contour  can  be  constructed  as  follows.  Let  F  be  the  median  triangle  in  the 
2-conduiL  Assigning  one  processor  to  each  triangle  in  the  /"-conduit  one  can  identify  the  set  S  of  all 
triangles  in  the  /* -conduit  which  interact  with  T.  This  set  is  an  interval  of  one  or  more  contiguous 
triangles  in  the  P-conduit;  let  Tj  and  T2  be  the  first  and  last  in  the  list  (these  may  be  the  same  trian- 
gle). Then  only  7,  and  the  triangles  above  it  in  the  /'-conduit  can  interact  with  triangles  above  T  in 
the  Q-conduit,  with  a  corresponding  observation  connecting  T2  with  the  triangles  below  T.  Thus  all 
pairs  of  interacting  triangles  can  be  identified  recursively  in  0(log(n))  parallel  steps.  Using  this 
information  it  is  straightforward  to  finish  constructing  the  contour  in  logarithmic  parallel  time. 

(F)  Finally  construct  the  cells  of  Vor{S)  by  fixing  up  the  cells  of  Vor (P),Vor  (Q)  as  appropriate. 
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5.  Intersections,  Partitions,  and  Related  Problems 

We  consider  the  following  problems:  (i)  detecting  whether  a  collection  of  line  segments  contains  any 
intersection,  (ii)  triangulating  a  simple  polygon,  (iii)  balanced  recursive  partition  of  a  polygon  (so  the  struc- 
ture of  the  partition  is  a  logarithmic  depth  binary  tree),  and  (iv)  computing  an  optimal  placement  of  watch- 
guards  in  an  art  gallery.  We  prove: 

Theorem  7.  The  above  problems  (i)-(iv)  are  all  in  NCiin  log(«)).  For  problems  (iii)  and  (iv),  if  we  begin 
with  a  triangulated  polygon  (or  art  gallery),  then  the  problem  in  in  NCiin). 

Note  that  using  straightforward  emulation  [Brent  (1974)],  the  NC\(n  log(/i))  algorithms  above  can  be  con- 
vened into  yVC  2(/i)-soludons.  These  time-bounds  have  recently  been  improved  to  0  ( log  (n))  [Atallah  and 
Goodrich  (1986),  Atallah,  Cole,  and  Goodrich  (1987)]. 

The  solutions  to  all  these  problems  share  a  common  basis,  which  is  the  construction  of  the  so-called 
(vertical)  trapezoidal  map  of  a  set  of  line  segments.  Let  V  be  a  collection  of  n  Unes  segments,  no  two  of 
which  intersect  except  perhaps  at  their  endpoints.  Following  [Lipski  and  Preparata  (1981)]  we  define  the 
trapezoidal  map  of  V  as  follows:  for  each  endpoint  p  of  a  segment  in  V,  let  Lp  denote  the  maximal  vertical 
line  segment  through  p  that  does  not  intersect  properly  any  segment  of  V.  (We  say  that  two  segments  inter- 
sect {Htjperly  if  their  relative  interiors  intersect)  As  is  apparent  from  the  illustration  in  Figure  5,  the  finite 
regions  of  the  trapezoidal  map  are  trapezoids. 


Figure  5.  Trapezoidal  map 

It  is  natural  to  combine  the  problem  of  segment  intersection  with  the  problem  of  computing  the  tra- 
pezoidal map:  we  will  give  an  algorithm  that  on  input  V,  detects  if  there  are  intersections  among  the  seg- 
ments of  V  and  if  not,  it  computes  the  trapezoidal  map  of  V.  For  brevity  in  exposition,  we  assume  no  two 
end  points  have  the  same  x-coordinate  and  in  particular  no  segment  of  V  is  vertical. 

Let  r  be  the  segment  tree  associated  with  the  projection  of  V  on  the  x-axis  [Bentley  (1977),  Bendey 
and  Wood  (1980)].  Recall  that  7  is  a  balanced  binary  tree  of  2n+\  leaves.  For  each  k  between  1  and 
2n+l,  the  Icth  leaf  of  T  from  the  left  is  associated  with  a  canonical  interval  on  the  x-axis:  this  interval  is 
delimited  by  the  k\h  smallest  x-coordinate  in  V  and  the  (;fc-(-l)st  smallest  one.   Similarly,  the  canonical 
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interval  of  every  internal  node  of  T  is  defined  recursively  as  the  union  of  its  two  children's  canonical  inter- 
vals. Let  Si s^  be  the  segments  of  V.  Each  node  v  of  the  segment  tree  is  assigned  a  set  of  indices, 

L(v),  called  its  node-set  as  follows:  let  /(j,)  be  the  projection  of  segment  s,  on  the  a:-axis.  We  say  that  a 
node  of  T  covers  I (Si)  if  its  canonical  interval  is  contained  in  /(^,)  but  the  canonical  interval  of  its  father  is 
not  contained  in  /(j,)  (or  else  the  node  has  no  father).  We  define  the  node-set  L  (v)  of  v  to  be  the  set  of  all 
indices  i  such  that  v  covers  I(Si).  It  is  straightforward  to  see  that  each  index  appears  in  <  21og(n)  nodes 
(at  most  two  nodes  per  level)  and  thus  that  the  storage  requirement  of  the  segment  tree  is  0  (n  log(n)). 

Let  us  return  to  our  original  problem.  Define  the  canonical  strip  of  a  node  v  of  7  to  be  the  infinite 
vertical  slab  formed  by  the  Cartesian  product  of  its  canonical  interval  and  (-«',-h>o).  The  intersection  of  the 
strip  with  the  segments  of  V  whose  indices  are  in  L(v)  forms  the  so-called  hammock  of  v,  denoted  H(v) 
[Chazelle  (1986)].  If  V  is  free  of  proper  intersections  then  the  hammock  consists  of  noncrossing  line  seg- 
ments joining  two  vertical  lines.  For  each  node  v,  define  the  set  W{v)  to  consist  of  segments  of  the  form 
s  r\Sy  where  5,  is  the  canonical  strip  of  v  and  s  intersects  the  interior  of  S,  and  s  has  at  least  one  endpoint 
in  the  closure  of  5,.  Again,  note  that  each  s  occurs  in  at  most  2  log  (n)  sets  W (v). 

First  we  consider  the  main  idea  of  the  algorithm.  The  problem  of  constructing  the  trapezoidal  map  of 

V  amounts  to  computing  the  segment  immediately  above  (below)  each  endpoint  in  V.  We  can  detect  if  the 
segments  in  H  (v)  intersect  by  sorting  the  segments  in  //  (v)  in  two  ways:  once  according  to  their  left  end- 
point  and  once  according  to  their  right  endpoint  (The  endpoints  lie,  respectively,  on  the  left  and  right 
boundaries  of  the  canonical  strip.)  If  these  sorting  orders  disagree,  we  at  once  detect  intersection.  Other- 
wise, for  each  segment  s  of  W(v),  we  can  insert  s  into  the  linearly  ordered  cells  of  the  canonical  strip  of  v 
in  two  ways:  by  inserting  s  according  to  its  two  endpoints.  Again,  s  e  W{v)  intersects  some  edge  in  //(v) 
if  and  only  if  the  two  cells  obtained  by  both  ways  of  inserting  s  into  the  canonical  strip  are  the  same.  If 
these  two  tests  for  H(v)  and  W{v)  pass  successfully  for  each  v,  then  we  can  conclude  that  the  segments  in 

V  do  not  intersect  (see  [Chazelle  (1986)]  for  a  prooO-  Note  that  it  is  insufficient  lo  just  do  these  tests  for 
//(v). 

Now  we  give  some  of  the  details  related  to  sorting.  The  tree  T  as  well  as  the  caironical  intervals  can 
be  computed  in  0  ( log  (n))  time,  using  sorting  on  n  processors.  Next  we  assign  one  processor  to  each  seg- 
ment of  V  and  compute  the  nodes  of  T  covering  each  segment  in  0  ( log  (n))  time  (using  binary  search  in 
the  tree).  At  this  point  we  can  assume  that  we  have  one  list  of  size  0  (n  log  (n)),  each  list  item  containing  a 
pair  of  the  form 

(z,v)  -  {segment  index,  covering  node). 
By  sorting  the  pairs  (/,v)  lexicographically  so  that  v  is  most  important,  we  get  the  pairs  having  a  common  v 
into  consecutive  locations  of  an  array.  We  may  regard  the  consecutive  locations  as  sublists  of  pairs 
corresponding  to  the  set  above  H  (v).  We  may  sort  the  edges  of  H  (v)  in  the  two  ways  described  above  and 
verify  that  no  intersection  occurs.  Next,  the  sets  W(v)  are  similarly  constructed  with  n  log(n)  pnxessors  in 
logarithmic  time.  We  may  then  insert  each  segment  s  o(W{v)  into  the  sorted  sublists  of  H  (v):  again  we  do 
this  in  two  ways  and  with  0  {n  log  (n))  processors  and  logarithmic  time.  So  now,  we  may  assume  that  we 
have  the  quadruples 

(p,V,Si,S2) 

where  p  is  an  endpoint  of  some  segment  s  inV,  W{v)  contains  a  portion  of  the  segment  s  which  terminates 
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in  p,  and  finally,  the  segment  ^i  (resp.  52)  come  from  the  hammock  edges  directly  above  (resp.  below)  p  in 
the  canonical  strip.  We  sort  again  to  collect  in  consecutive  locations  the  quadruples  involving  a  fixed  end- 
point  p.  By  computing  the  minimum  distance  among  the  Si  (resp.  ^2)  ^om  p,  we  obtain  the  edge  of  V 
directly  above  (resp.  below)  p.  This  completes  our  algorithm  for  the  trapezoidal  map.  To  prove  that  the 
algorithm  is  correct,  it  suffices  to  observe  that  by  finding  the  segment  immediately  above  (resp.  below)  an 
endpoint  p  will  necessarily  contribute  a  hammock  edge  which  lies  immediately  above  (or  below)  the  end- 
point  in  question. 

Note  that  the  data  structure  can  be  used  for  planar  point  location.  Better  yet,  it  allows  us  to  deter- 
mine the  segment  immediately  above  a  query  point  in  0  ( log  ^(/i))  time,  using  0  (n  log  (n))  space. 

In  the  next  two  subsections,  we  show  two  applications  of  the  trapezoidal  map. 

5.1.  Triangulating  a  Simple  Polygon 

To  triangulate  a  simple  n-gon  P,  we  first  compute  its  trapezoidal  map.  We  follow  [Chazelle  and 
Incerpi  (1984)]  in  the  following.  Our  goal  is  to  first  partition  P  into  monotone  polygcms.  To  do  this,  at 
each  reflex  comer  c  that  lies  in  the  relative  interior  of  a  vertical  edge  of  a  trapezoid,  we  introduce  a  diago- 
nal edge  from  c  to  the  comer  of  P  on  the  opposite  vertical  edge.  Such  diagonals  can  be  constructed  in  con- 
stant time.  This  now  partitions  P  into  a  collection  of  horizontally  monotone  polygons.  A  horizotally 
monotone  polygon  consists  of  an  upper  and  a  bwer  chain  of  edges,  which  meet  at  the  two  extremal  (left- 
most, rightmost)  comers.  A  chain  (upper  ot  lower)  is  trivial  if  it  consists  of  a  single  edge. 


Figure  6.  Horizontally  monotone  polygon  with  trivial  upper  chain 

We  show  how  to  partition  a  horizontally  monotone  polygon  Q  into  horizontally  monotone  subpo- 
lygons  with  a  trivial  upper  or  trivial  lower  chain,  as  follows.  Allocate  one  processor  to  each  edge  e  oi  Q. 
By  symmetry,  assume  that  e  forms  the  trivial  upper  chain.  If  the  region  of  the  plane  sufficiently  close  to  e 
and  below  e  is  outside  polygon  P  then  do  nothing.  Otherwise,  denoting  by  u^  and  ug  the  left  and  right  end- 
points  of  e,  we  (conceptually)  drop  the  vertical  line  segments  from  u^  and  u^  (resp.)  to  the  boundary  of  P. 
Let  these  Une  segments  meet  the  boundary  ot  P  at  v^  and  vg,  respectively.  (If  e  is  the  leftmost  (resp.  right- 
most) edge  of  its  upper  chain  then  u^  =  v£,  (resp.  ug  =  vg).)  Introduce  diagonals  from  u^  (ug)  to  the  left- 
most (rightmost)  vertex  of  P  between  v^^  and  Vg.  Compressing  the  vertices  of  each  monotone  polygon  into 
a  single  array  by  parallel  prefix,  this  produces  the  desired  partition  of  2  in  0  ( log  (/i))  time  using  q  proces- 
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sors,  where  q  denotes  the  number  of  vertices  of  Q  (Figure  6). 


So  now  we  may  assume  that  Q  =  [w^, . . .  ,w^]  is  a  horizontally  monotone  polygon  with  a  trivial 
upper  chain  w  i  w^,  and  lower  chain  given  by  the  polygonal  line  w  i , . . . ,  w^.  We  break  up  the  lower  chain 
into  "^  contiguous  segments  of  equal  size  and  compute  the  upper  hull  of  each  piece  (that  is,  the  upper 
chain  of  their  convex  hull).  This  is  similar  to  our  convex  hull  algorithm  in  section  3.  These  hulls  form  a 
continuous  path  that  partition  Q  into  a  connected  upper  portion  adjacent  to  w  i  w^  and  at  least  ^n  portions 
below  the  hulls.  Then  we  recursively  triangulate  the  parts  of  Q  below  each  upper  hull  (Figure  7). 


Figure  7.  Recursively  triangulate  below  each  of  the  "^  upper  hulls 

As  in  the  analysis  in  Section  3,  if  we  can  triangulate  the  part  of  Q  above  these  upper  hulls  in 
0  ( log^)  time  using  0  (q)  processors  then  the  complete  triangulation  of  Q  will  have  the  same  asymptotic 
complexity.  Once  again  by  assigning  one  processor  per  pair  of  upper  hulls  we  can  compute  all  pairwise 
tangents  in  logarithmic  time  (see  Section  3).  Next  we  build  up  a  partial  triangulation  of  the  region  of  Q 
between  these  upper  hulls  and  the  trivial  upper  chain:  1)  pair  up  consecutive  upper  hulls  (the  first  and 
second,  the  third  and  fifth,  etc.)  and  for  each  pair  determine  the  unique  tangent  between  the  two  upper 
hulls,  adding  it  into  the  triangulation;  2)  iterate  on  this  process  until  there  is  only  one  upper  hull  lefL  Step 
(1)  can  be  implemented  in  constant  time:  lo  see  this,  assign  one  processor  per  vertex,  per  ciurent  upper  hull 
and  per  precompuied  tangent  At  the  beginning  of  each  new  iteration  we  assume  that  the  processors 
assigned  to  vertices  which  are  no  longer  on  the  current  upper  hulls  are  de-activated.  Then  in  0  (1)  time  the 
vertices  of  the  current  upper  hulls  are  each  labelled  with  a  i.d.  number  unique  to  the  hull.  In  this  way  each 
(tangent)  processor  can  then  decide  immediately  if  it  is  relevant,  that  is,  if  it  connects  two  hulls  in  a  given 
pair.  If  so,  then  the  tangent  is  added  into  the  triangulation  and  the  processor  of  every  vertex  lying  strictly 
below  the  tangent  is  de-activated.  Tangents  which  crosses  the  newly  added  edges  of  the  triagulation  also 
have  their  processors  de-activated.  The  iteration  can  then  proceed.  When  the  iterations  termini,^  we  have 
the  situation  of  figure  8. 
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Figure  8.  Partial  triangulalion  of  region  above  the  upper  hulls 

It  is  easy  to  see  that  the  partial  triangulation  leaves  behind  polygonal  regions  with  very  simple 
characteristics:  such  polygons  2'  =  (u  i , .  . . ,  u*]  are  horizontally  monotone  with  a  trivial  upper  chain  m  i  u^, 
and  there  is  a  unique  lowest  vertex  ui(\  <l  <k)  such  that  the  y-coordinate  of  u,  is  decreasing  as  i  increases 
from  1  to  /,  and  the  increasing  as  i  increases  from  /  to  k.  Furthermore,  the  subpolygons  {mj,  . . .  ,u,)  and 
{u,, . . . ,  ut )  are  also  horizontally  monotone  with  a  trivial  lower  chains  u i  ui  and  uiu^,  respectively  (Figure 
9).  Since  they  are  also  convex,  it  is  now  an  easy  exercise  to  triangulate  each  Q'  using  k  processors  and 
0  ( log  (k))  time. 


Figure  9.  Special  monotonic  polygons  after  partial  triangulation 

In  conclusion,  we  have  given  a  algorithm  ior  triangulating  a  simple  n-gon  in  0  ( log  (n))  time,  using  n 
processors  if  the  vertical  trapezoidal  map  is  available,  or  using  0  (n  log  (n))  processors  otherwise. 

5^.  Poiygon-Cutting,  Recursive  Decompositions  and  Applications 

The  polygon-cutting  thecrem  states  that  any  simple  polygon  has  a  diagonal  which  cms  it  into  two 
roughly  equal  pieces.  In  its  simplest  form,  it  says  that  an  n-gon  can  always  be  cut  up  into  polygons  of  at 
most  2(/i/3)-(-l  vertices.  As  was  shown  in  [Chazelle  (1982)]  this  property,  though  quite  simple,  is  of  great 
importance,  because  like  other  separator  theorems  it  sets  the  stage  for  divide-and-conquer  algorithms.  We 
will  illusffate  this  with  the  problem  of  placing  guards  in  an  art  gallery.  First  we  review  how  the  polygon- 
cutting  theorem  is  used. 

In  many  problems  on  simple  polygons,  it  is  useful  to  divide  a  polygon  roughly  evenly  for  separate 
processing,  and  to  recursively  proceed.  This  amounts  to  constructing  a  'decomposition  tree'  that  is  bal- 
anced (i.e.  logarithmic  depth),  where  each  node  v  stores  a  diagonal  of  the  polygon  (representing  a  partition) 
and  the  left  and  right  subtrees  below  v  represent  the  two  polygons  defined  by  this  partition;  the  leaves 
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represent  triangles.  Note  that  the  set  of  all  diagonals  stored  in  the  decomposition  tree  forms  a  triangulation 
of  the  polygon.  The  polygon-cutting  theorem  implies  the  converse:  every  triangulation  can  be  organized 
into  a  balanced  decomposition  tree.  The  graph-theoretic  dual  of  the  triangulation  is  a  tree  whose  nodes  and 
edges  are  in  bijective  correspondence  with  the  triangles  and  the  diagonals,  respectively.  If  a  diagonal  parti- 
tions the  boundary  of  P  into  two  polygonal  Unes  of  respectively  p  and  n-p  edges  then  removing  the 
corresponding  edge  from  the  tree  will  leave  two  trees  of  respectively  p-1  and  n-p-\  nodes.  What  is  now 
needed  is  an  efficient  parallel  algorithm  for  computing  the  size  of  the  subtrees  adjacent  lo  every  edge  in  the 
tree.  This  can  be  reduced  lo  computing  the  number  of  descendents  of  every  node  in  a  rooted  tree,  which 
we  have  already  used  in  our  algorithm  for  Voronoi  diagrams  (section  4).  From  this  we  conclude  that  bal- 
anced decomposition  trees  can  be  implemented  in  NC\{n)  if  the  polygon  is  already  triangulated,  and  in 
NC\{n  log  (n))  otherwise. 

In  [Chvdtal  (1975)]  it  is  shown  that  as  long  as  P  is  simple  it  is  always  possible  to  place  [n/3|  watch- 
gxiards  inside  P  to  keep  the  whole  polygon  ('art  gallery')  in  check.  This  means  that  with  Ynll^  carefully 
chosen  points  every  point  in  f  is  visible  from  at  least  one  of  the  chosen  points.  A  simple  proof  follows 
from  the  property  that  any  triangulation  of  f  is  3-colorable.  Being  a  tree,  the  dual  graph  of  a  triangulation 
of  P  has  al  least  one  leaf.  This  means  that  if  P  has  more  than  three  vertices,  its  triangulation  contains  at 
least  one  triangle  abc  adjacent  to  the  rest  of  P  through  ab  only.  Removing  this  triangle  allows  us  to  3-color 
the  remainder  recursively,  which  then  leaves  one  color  available  for  vertex  c.  A  solution  to  the  art-gallery 
problem  then  results  from  placing  a  guard  at  each  vertex  colored  with  the  least-used  color. 

To  overcome  the  inherendy  sequential  nature  of  this  algorithm,  we  first  assume  that  a  decomposition 
tree  T  for  the  polygon  is  available.  Our  goal  is  to  3-color  a  simple  polygon  P.  A  3-coloring  of  P  is 
represented  by  storing  at  each  edge  of  T  a  permutation  of  the  colors  c i ,  ci,  C3,  and  at  the  leaves  of  T,  an 
arbitrary  3-coloring  of  the  vertices  of  the  triangle.  Note  that  a  vertex  of  P  may  appear  in  more  than  one  tri- 
angle but  we  make  no  assumption  that  the  color  assigned  to  its  various  appearances  at  the  leaves  are  con- 
sistent The  consistency  question  is  resolved  by  the  permutations  at  the  edges.  Let  us  describe  how  the 
color  of  any  vertex  v  of  /*  is  obtained  in  this  coding  scheme.  Say  v  appears  in  the  triangle  at  a  leaf  /  of  the 
decomposition  tree  T.  If  v  is  'locally'  assigned  a  color  c  at  /,  then  the  color  of  v  in  the  'global  coloring'  is 
obtained  by  applying  the  permutations  found  along  the  path  from  /  to  the  root  of  T.  We  must  ensure  con- 
sistency, i.e.,  different  leaves  containing  v  will  lead  to  the  same  global  color  for  v.  This  is  simple.  Let  us 
say  that  a  node  u  of  T  is  globally  consistent  if  for  all  vertices  v  of  P,  whenever  two  leaves  below  u  assigns 
local  colors  to  v,  then  the  path  from  these  two  leaves  to  u  will  assign  the  same  global  color  to  v.  Assume 
that  all  nodes  at  distance  greater  than  k  from  the  root  are  globally  consistent  Let  u  be  a  node  at  distance  k 
from  the  root,  and  u^,  ug  be  the  two  children  of  u.  It  is  easy  to  assign  permutations  to  the  two  edges  from  u 
to  its  children  so  that  u  is  globally  consistent  (To  see  this,  not  that  the  possible  cause  for  inconsistency  is 
at  the  two  vertices  incident  to  the  diagonal  of  /*  at  u.)  Starting  firom  the  nodes  furthest  away  fronr  the  root, 
in  0  ( log  (n))  steps,  we  can  compute  the  permutations.  In  another  O  ( log  (/i)  steps,  with  one  processor  per 
leaf  of  T,  we  can  assign  the  global  color  to  each  v  in  P. 

6.  Polygon  Optimization  using  a  Back-and-Forth  technique 

There  has  been  considerable  interest  recently  in  optimization  problems  where  the  solution  space  is  a 
suitable  class  of  polygons  (see  [Chang  (1986)1).   One  such  class  of  problems  is  the  following,  for  fixed 
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k>3.  Given  a  convex  n-gon  P,  determine  the  minimum-area  (resp.  maximum-area)  ifc-gon  that  cir- 
cumscribes (resp.  is  inscribed  in)  P.  First  consider  the  circumscribing  problem.  In  the  sequential  setting, 
[O'Rourke  et  al  (1985)]  give  a  linear  time  algorithm  for  yfc  =  3  while  [Aggarwai.Chang  &  Yap  (1985), 
Chang(1986)]  give  an  C>(n^log(/i)log()t))  algorithm  for  k>4.  Recently,  a  factor  of  login)  has  been 
shaved  off  this  result  using  a  general  technique  in  [Aggarwal,  Klawe,  Moran,  Shor  &  Wilber  (1986)]. 
Since  the  solution  for  i  >  4  is  at  least  quadratic  time,  we  cannot  have  an  NCX(n)  solution  without  improv- 
ing the  sequential  result  Thus  we  focus  on  the  minimum  circumscribing  triangle  problem.  Our  main  result 
is  given  next  the  stated  nmtime  is  achieved  using  what  we  have  dubbed  'back-and-fcMth  subdivision,'  a 
technique  new  in  computational  geometry. 

Theorem  8.  The  problem  of  computing  a  minimum  circumscribing  triangle  for  a  convex  n-gon  is  in 
NC\(n). 

For  convenience  we  refer  to  the  n  edges  of  a  given  convex  polygon  P  by  their  indices  in  clockwise  order 
firom  0  to  n-1.  Also,  given  an  edge  i,  we  denote  by  a(/)  the  last  edge  in  clockwise  order  from  i  which 
forms  a  positive  angle  with  the  line  through  i:  put  another  way,  a(/)  is  the  unique  edge  such  that  its  second 
terminal  (in  clockwise  order)  supports  the  tangent  to  P  parallel  to  /  and  opposite  /  (there  is  an  easy  special 
case  when  this  tangent  lies  along  an  edge).  We  exhibit  a  ^n  divide  and  conquer  strategy  (cf.  Section  3) 
which  relies  on  the  following  geometric  characterization  developed  in  [Klee  and  Laskowski  (1985), 
O'Rourke  et  al  (1985)]. 


Figure  10.  The  edges  {/,/-i-l, . . .  ,a(j))  o(P  and  T{i,j). 

We  use  the  following  known  properties: 

(i)  Among  the  minimum  circumscribing  triangles  of  P,  at  least  one  such  triangle  has  two  of  its  sides 
'flush'  (where  a  side  is  flush  if  it  lies  along  an  edge  of  P). 

(ii)     The  midpoint  of  each  side  of  any  minimum  circumscribing  triangle  of  P  must  touch  P. 

(iii)  Given  i  and  y,  where  i<j<a{i),  let  7"(j,  j)  denote  the  triangle  with  two  sides  flush  with  i  and  j  such 
that  the  third  side  has  its  midpoint  q  touching  P  and  such  that  z,  ;  and  q  occur  in  clockwise  order  on 
the  boundary  of  P.  We  write  m (/,_/)  for  the  edge  of  P  that  touches  q.  In  general  the  third  side  of 
T{iJ)  is  not  flush  so  there  will  be  two  choices  for  m(ij);  we  will  conventionally  choose  the  larger 
index  (mod  n).  For  each  /,  there  are  at  most  two  indices  j*(i)  such  that  the  midpoint  of  the  second 
side  of  T{i,  f*  (/))  (i.e.,  the  side  that  is  flush  with  _/'*  (/))  actually  touches  P.  If  there  are  two  choices 
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for  y*  {/),  we  choose  the  edge  that  has  the  larger  index  (mod  n).  Let  m*  (/)  refer  to  the  edge 
miij*  (/)),  and  T  (i)  refer  to  7"(/,  ;'*  (/)).  We  say  edge  j*{i)  (resp.  m*  (0)  is  left  critical  (resp.  right 
critical)  with  respect  to  i.  We  will  call  triangle  T*  (i)  a  locally  minimum  triangle  with  respect  to  edge 
i.  Note  that  the  midpoint  of  the  edge  of  T*  (/')  flush  with  i  need  not  touch  P. 

(iv)  Monoionicity  property.  An  edge  ye/'  (i  *j)\s  called  low  {high,  resp.)  with  respect  to  /  if;  precedes 
y*(j)  0  succeeds  j*{i),  resp.)  in  a  clockwise  traversal  of  /'.  By  definition  then,  for  each  /,  the 
sequence  of  edges  (/,i+l, . . .  ,ot(i))  consists  of  a  sequence  of  low  edges,  followed  by  a  left  critical 
edge,  followed  by  a  sequence  high  edges.  The  monotonicity  property  is  as  foUows.  Suppose  thatp 
is  the  midpoint  of  the  side  of  T(iJ)  lying  along  j  and  p  is  not  on  ;:  if  p  and  the  side  containing  i  are 
on  the  same  side  of  j  then  ;  is  high,  and  if  they  are  on  opposite  sides  then  j  is  low,  with  respect  to  /. 
(Note  that  this  deduction  is  made  without  knowledge  of  the  left  or  right  critical  edges  of  i.) 

(v)      Interspersing  property.  If  /'  <  i'  then  j*  (/)  S  j*  (0  and  m*  (/)  <,  m*  (0- 

Property  (v)  is  a  weakened  form  of  that  derived  by  O'Rourke  et  al.  (1986)  but  is  sufficient  for  our 
purposes.  Given  i  and  j  one  can  compute  m*  (ij)  using  binary  search  in  logarithmic  serial  time;  thus,  given 
i,  the  monotonicity  property  lets  us  use  binary  search  to  compute  j'(i)  in  0  (log^(n))  serial  time.  Thus  n 
processors  can  compute  the  locally  minimum  triangles  T*{i)  for  each  i.  By  property  (i),  the  globally 
minimum  triangle  will  be  found  among  these  T*{i).  This  can  be  found  in  0(log(n))  further  steps.  We 
shall  see  that  this  can  be  improved  to  yield  an  0(log(n))  parallel  algorithm  using  what  we  shall  call 
'back-and-forth  subdivision.' 

We  define  partial  right  inverses  for  the  two  functions  j*  and  m*:  _/"'  (k)  is  the  unique  i  such  that 
j*{i-l)<k<j*{i)  (in  clockwise  order);  m~\k)  is  defined  the  same  way.  (The  functions  j*  and  m*  are  not 
usually  bijective,  but  if ;(:  is  in  the  range  of  j*  then  y*0~'  ik))=k\  similarly  for  m~' .)  The  monotonicity  pro- 
perty (iv)  shows  that  these  quantities  are  well-defined  and  furthermore,  given  an  interval  /  of  length  s  con- 
taining j~^(k),  if  /*(/)  is  known  for  all  i  in  that  interval,  then  a  single  processor  can  compute  j~^(.k)  in 
0  ( log  (s))  time.  The  idea  is  that  we  shall  compute  approximate  values  for  functions  f*  and  m*  and  for 
their  inverse  functions  ;"'  and  m~' .  We  can  refine  an  estimate  oij*{i)  by  using  information  about  j~^  and 
m*  and  performing  binary  search  in  the  interval  known  to  contain  the  true  value.  The  other  approximations 
can  be  refined  simultaneously.  Repeating  this  process  sufficiently  often,  namely,  log  log(n)  times,  we 
compute  all  these  functions  exactly.  A  single  phase  of  the  process  involves  knowing  for  each  k  an  interval 
a  ■  ■  ■  b  of  length  s,  say,  such  that  a<j~^(k)<b.  In  this  case  we  say  that  the  interval  [a,b]  brackets  j*{k). 
Using  sufficiently  many  processors  (n  overall),  we  can  improve  this  estimate  to  an  interval  of  length  "^s: 
the  same  goes  for  calculating  estimates  of  the  other  functions.  The  phase  takes  0  ( log  (s))  parallel  steps. 

For  simplicity  we  shall  assume  that  /i  is  a  power  of  two.  We  initially  choose  5  as  the  maximum 
integer  <a  of  the  form 

where  t  is  an  integer,  (specifically,  i^  log  2(  log  2(1))!  ).  This  avoids  rounding  difficulties.  This  initial  value 
of  5  lies  between  "*/«  and  n. 

Lemma  9.  Let  s  be  a  perfect  square  perfectly  dividing  n  (which  is  assumed  to  be  a  power  of  2).   For 
0<i<n/^s-l  let  s,  be  the  i^s-th  edge  in  clockwise  order  around  the  polygon  P.  Suppose  that  to  each  k  is 
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associated  four  intervals  J*{k).  J~^{k),  M*{k).  and  M~^{k)  respectively  bracketing  j*(k),  j~^(k),  m*{k), 
and  m~^(k),  where  each  interval  is  of  the  form  [Si-JJ,s^i+i)-^].  In  other  words,  the  given  bracketing  inter- 
vals are  of  length  s+l  and  their  boundaries  are  multiples  ofs.  Then  with  n  processors  all  of  these  intervals 
can  be  reduced  to  intervals  of  the  form  [^i.^i+J.  in  0{log(s))  parallel  steps. 

Proof,  For  0<i<n/^  we  assign  \y  processors  to  compute  in  0(\og(s))  parallel  steps  the  correct 
subinterval  [j,,j,+i]  of  7*(j,)  bracketing  j*(s,).  The  f-th  processor  uses  binary  search  on  A/*(j,)  to  deter- 
niine  whether  s,<j*(Si),  i.e.,  to  determine  whether  s,  is  low,  critical,  or  high  with  respect  to  Si.  Thus  the 
correct  interval  bracketing  j*  {s,)  is  located.  Reassigning  the  given  processors  to  this  interval  7* (i,)  is  com- 
puted exactly.  This  exact  information  about  f*{Si)  at  these  n/ ^5  sample  sides  j,  enables  us  to  bracket 
y"'  (k)  for  all  sides  kofP,  since 

Si<r\k)<Si^i 
if  and  only  if 

j*isi)<k<r(s,^,), 

so  this  can  be  checked  in  time  0  ( log  (s))  with  one  processor  reassigned  to  each  edge  k,  using  binary  search 
on  the  subsequence  of  edges  5,  contained  in  the  interval  J"'  (;k)  0*  (s,)  is  known  exactly  for  all  i). 

Thus  we  have  j*(,Si)  computed  for  all  i  and  y~'  (k)  correctly  bracketed  for  all  k.  Next  for  each  sample 
side  Sk  we  assign  ^s  processors:  one  is  assigned  to  each  side  /  in  the  interval  /"'  (;fc)  (which  is  now  of  the 
form  [5,,  J(+i  ]).  The  processor  assigned  to  /  queries  M*  (i)  and  determines  whether  k  is  low,  critical,  or  high 
with  respect  to  i.  By  the  monotonicity  property  this  enables  j~^  {st)  to  be  determined  exactly.  When  this 
computation  is  finished,  /*(;')  can  be  recomputed  for  every  edge  i  just  as  /~'(0  was  computed.  The 
refinement  of  our  estimates  of  M*  (k)  and  Af "'  (k)  is  achieved  by  symmetric  methods.  Q£.D. 

The  initial  conditions  for  applying  the  lemma  are  slightly  different  since  in  general  n>s.  The  aim 
then  is  to  achieve  initial  bracketing  intervals  of  length  r.  however,  essentially  the  same  ideas  apply. 
Repeating  the  refinement  phase  log  log(/i)  times  all  quantities  are  computed  exactly.  Notice  that  since  all 
the  tasks  assigned  involve  easily-computed  subsequences  and  subintervals  of  0  ■  •  •  n-1,  processor  assign- 
ment presents  no  difficulty  in  this  algorithm.  It  is  clear  that  this  method  takes  overall  parallel  time 
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SO  ( log  (n^  ))  =  YP  (^^^^)  =  O  ( log (n)). 
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Additional  Remarks.  Boyce  et  al  (1985)  have  shown  that  some  of  the  properties  similar  to  those 
given  for  the  minimum  circumscribing  triangle  also  hold  for  the  maximum  inscribed  triangle.  It  is  not  clear 
if  the  back-and- forth  technique  can  be  extended  here  since  the  monotonicity  property  (iv)  is  not  available. 
However  an  (9  ( log  ^(n))  solution  is  possible.  For  the  dual  problem  of  computing  the  maximum  inscribed 
triangle,  simple  divide  and  conquer  gives  a  solution  using  (9(log^(n))  time  and  n  processors.  In  fact, 
unlike  the  minimum  circumscribed  /t-gons  (for  k  >  4),  such  a  divide  and  conquer  technique  yields  an 
0{k\og^{n))  time  for  computing  a  maximum  inscribed  ;k-gon  using  n  processors.  Finally,  Toussaint 
(1983)  has  shown  that  the  minimum  circumscribing  rectangle  for  a  given  convex  n-gon  can  be  found  in 
0{n)  sequential  dme.  Using  the  geometric  characterization  given  by  Freeman  and  Shapira  (1975),  it  is 
easy  to  compute  such  a  rectangle  in  0  ( log  (/i))  time  using  only  nl  log  (n)  processors. 

In  fact,  the  diameter  of  a  convex  polygon  can  also  be  computed  in  0  ( log  (n))  time  using  n  I  log  (n) 
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processors,  as  the  following  brief  outline  suggests.  If  we  can  compute  for  each  i  its  antipodal  vertex  5(i) 
(the  endpoint  of  (x(/))  it  is  straightforward  to  compute  the  distance  of  5(/)  from  the  line  through  /  for  each  i, 
and  then  compute  the  minimum  of  all  these  distances  with  the  stated  resources.  We  take  n/  log  (n)  proces- 
sors and  assign  each  processor  to  a  group  of  log(n)  successive  sides  of  the  polygon.  Each  processor  first 
computes  the  antipodal  vertex  for  the  first  side  in  its  alotted  group  (this  involves  unimodular  search  in  loga- 
rithmic serial  time),  and  then  begins  to  scan  in  a  clockwise  direction  along  its  group  of  sides.  Its  'primary 
task'  is  to  compute  5(/)  for  all  sides  within  this  group.  Going  from  side  i  to  j+l,  the  corresponding  antipo- 
dal vertex  may  either  remain  the  same,  increase  by  1 ,  or  iiKrease  by  mwe  than  1 .  In  the  latter  case  the  gap 
between  5(0  aixl  5(/+l)  could  be  considerably  more  than  log(/i),  and  the  processor  might  not  achieve  its 
primary  task.  However,  in  this  case,  if  k  denotes  the  side  beginning  at  5(0,  we  know  that  5{k)  is  the  vertex 
ending  the  j'-th  edge.  Thus  the  processor  can  begin  a  'secondary  task,"  namely  computing  50)  for 
j=k,k+l,  •  ■  ■  until  it  either  reaches  the  end  of  the  group  containing  k  or  5(J)  increases,  in  which  case  it  can 
determine  5(i+l)  and  revert  to  its  primary  task.  It  is  easy  to  show  that,  for  every  i,  5(i)  is  computed  by 
some  processor,  either  as  a  primary  or  secondary  task. 


7      Convex  hull  in  3  dimensions  and  related  structures. 

In  this  section  we  shall  discuss  the  problem  of  computing  the  convex  hull  H(S)  of  a  set  S  of 
points  in  three  dimensional  space.  The  method  originally  described  in  the  FOCS  paper  ran  in 
parallel  time  0(log''(n)),  but  the  method  to  be  presented  here  is  in  NC^in).  Thus  we  now  match 
the  asymptotic  efficiency  of  Anita  Chow's  method  [1980],  and  we  do  not  need  an  optimaJ  parallel 
sorting  algorithm.  As  with  the  2-dimensional  Voronoi  diagram  problem  the  algorithm  is  based 
on  a  serial  divide-and-conquer  technique  alreauiy  known  in  the  literature.  We  shall  employ  a 
new  data-structure,  which  has  independent  interest,  to  solve  queries  of  the  foUowing  kind:  given 
a  line  L  and  a  polytope  K,  to  determine  whether  L  intersects  K  and  if  not  to  return  a  plane 
through  L  tangent  to  K .  Such  queries  will  be  called  line-queries;  they  will  be  discussed  in  detail 
later.  Let  us  state  the  main  result  of  this  section: 

Theorem  10  There  is  an  NC^in)  algorithm  for  constructing  the  convex  hull  H{S)  of  a  set  S 
of  n  points  in  3-space,  which  only  uses  sorting  at  a  preprocessing  stage  and  hence  can  use  any 
sorting  algorithm  tn  NC^{n)  or  better.  The  algorithm  produces  as  a  by-product  a  data- structure 
suitable  for  answering  line-queries  about  H{S)  in  serial  time  0(log  (n)). 

Recall  the  serial  divide-and-conquer  algorithm  [Preparata  and  Hong  (1977)]  to  compute  the 
convex  huU  of  S:  the  set  of  points  S  is  evenly  split  into  two  sets  by  a  horizontal  plane  R,  which 
does  not  contain  any  point  in  5;  Let  P  and  Q  denote  the  points  in  5  above  and  below  R, 
respectively.  Recursively  compute  H{P)  and  H{Q);  compute  the  sleeve  joining  the  two  hulls, 
i.e.,  the  chain  of  faces  tangent  to  the  two  hulls  and  meeting  three  or  more  points  of  S;  then 
merge  H{P)  with  H{Q)  along  the  sleeve  to  form  the  convex  hull  H{S). 

For  the  rest  of  this  section  R  wOl  continue  to  denote  a  horizontal  plane,  which  we  shall  call 
the  separating  plane,  partitioning  S  into  upper  and  lower  parts  P  and  Q  respectively.  Consider 
that  part  of  the  sleeve  on  and  above  the  separating  plane  R.  One  might  imagine  that  it  was 
homeomorphic  to  a  cylinder  with  boundary,  but  sleeve  faces  can  meet  one  another  above  the 
separating  plane,  so  it  can  have  several  holes.  Each  partial  sleeve  face  has  an  edge  contained  in  R, 
and  two  edges  Incident  on  this  edge;  all  other  edges  on  this  partial  face  will  be  called  seam  edges 
and  given  a  left-to-right  (i.e.,  clockwise)  orientation.  These  chains  of  seam  edges  join  together 
to  form  what  will  be  called  the  upper  seam  of  the  H{S).  Thus  the  upper  seam  has  the  form  of 
a  directed  Eulerian  circuit,  but  vertices  may  be  visited  more  than  once  and  the  same  edge  of 
H{S)  may  occur  with  both  orientations  along  the  seam."*  The  same  edge  e  in  H{P),  say,  will  be 
visited  in  both  directions  by  the  seam  if  both  the  front  and  bcick  tangent  planes  (see  Definition 
11)  through  that  edge  to  H{Q)  are  tangent  to  H{P);  it  is  easy  to  construct  examples  where  this 
happens.  The  lower  seam,  occurring  below  the  separating  plane  R,  is  defined  similarly. 

Our  algorithm  is  simple  in  principle  though  the  development  is  lengthy:  it  depends  on  com- 
puting the  seams  efficiently  at  each  stage.  To  compute  the  upper  seam,  say,  cissign  one  processor 
to  each  edge  of  H{P),  to  peform  what  will  be  called  a  'Une  query,'  relative  to  H{Q).  Consider 
a  fixed  edge  of  H{P);  let  L  be  the  Une  containing  it.  If  the  edge  belongs  to  the  upper  seam  it 
is  necessary  (a)  that  L  not  meet  the  interior  of  H{Q),  and  (b)  one  (or  conceivably  both)  of  the 
two  planes  through  L  tangent  to  H{Q)  should  not  intersect  the  interior  oi  H{P).  To  determine 
whether  a  given  plane  passing  through  L  \s  tangent  to  H{P)  can  be  accomplished  in  constant 
time  by  inspecting  the  two  faces  of  H{P)  meeting  L. 

Once  the  edges  in  the  seam  have  been  calculated,  their  cyclic  connection  order  can  be  as- 
certained in  logarithmic  time  by  a  list-ranking  process.  The  full  structure  of  the  sleeve  can  be 
deduced  once  both  seams  have  been  constructed  by  implementing  what  is  essentially  a  merging 
process,  easUy  accompUshed  in  logarithmic  time.  Therefore  we  shall  now  consider  these  line 
queries  more  closely. 

Let  A"  be  a  convex  polytope  and  L  an  oriented  hne  not  meeting  the  interior  oi  K .  There  are 
two  half-planes  through  L  tangent  to  K  (they  bound  the  convex  hull  of  LuK:  let  us  ignore  the 

*The  possibility  that  the  seeun  can  have  self-intersections  complicates  both  the  sequential  algorithm  and  ours. 
We  are  indebted  to  Herbert  Edekbrunner  (private  conununication)  for  drawing  our  attention  to  this  possibility. 
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Figure  11:  Illustrating  the  front  and  back  tangent  planes 

trivial  case  where  K  and  L  are  contained  in  the  same  plane).  It  is  helpful  to  distinguish  these 
planes  in  a  canonical  way. 

Definition  11  With  K  and  L  as  above,  consider  the  set  of  all  half-planes  bounded  by  L  which 
intersect  K.  All  such  half-planes  can  be  parametrized  by  a  rotational  angle  9  about  L.  Assmntng 
thai  L  does  not  intersect  (the  interior  of)  K,  these  angles  9  sweep  out  an  interval  [^o.^i]  of 
angles  where  we  can  distinguish  these  two  bounding  angles  (since  L  is  oriented)  by  requiring  thai 
6i  >  9q  in  the  anticlockwise  sense,  and  9i—do  <  tt.  These  bounding  angles  define  the  two  planes 
through  L  tangent  to  K,  and  we  call  the  one  at  angle  9q  ike  back  tangent-plane  to  K  through 
L,  and  the  other  the  front  tangent-plane.  The  line-query  for  L  and  K  is  to  determine  whether 
L  meets  K  transver sally,  and,  if  not,  to  compute  the  front  and  back  tangent  planes  through  L  to 
K. 

If  we  visuahze  L  and  A'  so  that  L  is  in  the  plane  of  the  page,  and  oriented  upwards,  and  K 
is  to  its  left,  then  the  front  and  back  planes  will  be  in  front  of  and  behind  K  respectively  (see 
figure  11). 

Setting  up  the  means  to  solve  such  queries  will  occupy  most  of  this  section.  To  begin  with, 
let  us  consider  a  rather  simple  query  which  we  call  a  point-query. 

Given  a  plane  T  meeting  a  corner  p  of  a  convex  polyhedron  A',  to  determine  whether 
T  is  tangent  to  K ,  and  if  not  to  give  the  two  faces  incident  to  p  which  T  cuts. 

This  is  solved  by  providing  a  test  polygon  at  p,  namely  a  polygon  which  forms  the  boundary  of 
the  region  of  intersection  of  a  fixed  plane  M  separating  p  from  all  the  other  corners  of  A'.  Clearly 
the  plane  T  cuts  K  if  and  only  if  it  cuts  this  test  polygon,  and  if  so  then  the  edges  where  it  cuts 
the  polygon  identify  the  faces  incident  to  p  which  it  cuts.  The  corners  of  the  test  polygon  are 
unimodal  with  respect  to  height  above  (or  below)  the  plane  T,  so  it  is  easy  then  to  solve  the 
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queries  in  logarithmic  serial  time  by  a  variant  of  binary  search.  The  problem,  therefore,  is  to 
construct  a  test  polygon  in  parallel,  with  one  processor,  say,  assigned  to  each  edge  of  K  incident 
to  p.  Indeed,  it  is  eeisy  to  construct  such  a  polygon  once  at  least  one  line  A''  has  been  found 
which  meets  the  interior  of  K  at  p,  since  then  a  suitable  separating  plane  normal  to  this  line 
can  be  found  (take  all  the  projections  of  the  incident  edges  on  N ,  and  find  the  one  closest  to  p 
in  logarithmic  time;  let  r  be  its  distance  from  p;  then  choose  the  plane  passing  perpendicularly 
through  A^  at  half  this  distance).  To  find  a  line  entering  K  at  p,  simply  choose  three  corners 
adjacent  to  p  and  let  q  be  their  average;  then  pq  defines  a  suitable  line.  Summarizing, 

Lemma  12  Given  a  convex  polyhedron  K  with  I  comers,  it  is  possible  to  construct  test  polygons 
for  each  comer  in  0{\og{£))  parallel  time  overall  with  t  processors. 

The  following  lemma  shows  where  point-queries  could  arise  in  constructing  the  convex  hull. 
Its  simple  proof  is  omitted. 

Lemma  13  Let  T  be  a  plane  which  meets  the  sleeve  of  H{S)  at  some  comer  p  but  does  not 
intersect  any  face  of  the  sleeve.  Then  T  can  meet  the  interior  of  H{S)  if  and  only  if  it  meets  the 
interior  at  p,  and  in  this  case  the  polygonal  region  of  intersection  of  T  with  H{S)  is  contained 
in  H{P)  if  p  is  in  the  upper  seam  and  in  H(Q)  if  p  is  in  the  lower  seam. 

The  rest  of  this  section  is  organized  as  follows.  We  introduce  the  notion  of  'seam  polytope,' 
and  of  locahzing  a  query  to  a  siib-polytope.  We  study  the  sequencing  and  nesting  properties  of 
the  seam  on  a  seam  polytope,  and  use  them  to  give  that  part  of  the  polytope  above  the  seam 
the  structure  of  a  forest  of  rooted  trees.  We  introduce  the  idea  of  a  'near  split'  on  a  polytope 
which  can  be  used  to  localize  a  line-query  on  the  polytope.  Next  we  define  the  'core  and  lobes' 
of  a  seam  polytope,  and  show  how  to  localize  a  line-query  to  one  of  these  lobes.  Next  we  show 
how  to  decompose  a  lobe  recursively  by  a  sequence  of  near-spUts,  in  a  balanced  manner,  thereby 
completing  the  construction  for  solving  hne-queries.  Finally  we  show,  given  a  description  of  a 
sleeve  with  its  incident  edges  and  faces  (in  cyclic  order),  how  to  build  the  two  associated  seam 
polytopes  (a  seam  polytope  is  obtained,  roughly  speaking,  by  discarding  irrelevant  faces  from 
H{S):  the  resulting  polytope  will  only  resemble  H{S)  close  to  the  seam,  so  it  is  necessary  to 
build  it  from  scratch). 

7.1      Definition  and  cliaracteristics  of  a  seam  polytope 

For  each  of  the  two  seams  of  the  sleeve  as  considered  above  we  shall  construct  a  'seam  polytope,' 
a  bounded  convex  polytope  which  has  a  distinguished  circuit  of  directed  edges,  called  the  seam, 
matching  the  corresponding  seam  of  the  sleeve. 

Definition  14  A  seam  polytope  is  a  convex  3-dimensional  polytope  with  two  distinguished  hor- 
izontal faces,  the  base  and  the  top,  and  the  other  faces  partitioned  into  two  sets,  called,  for 
mnemonic  purposes,  the  green  and  blue  faces  respectively.  The  edges  along  the  top  and  base 
are  called  top  and  base  edges  respectively.  No  green  face  meets  the  base  and  no  blue  face  meets 
the  top.  Each  blue  face  meets  the  base  in  a  proper  edge  and  is  either  triangular  or  trapezoidal 
in  shape.  The  seam  vertices  and  edges  are  those  vertices  and  edges  bounding  blue  faces  but  not 
meeting  the  base.  Edges  connecting  the  seam  to  the  base  are  called  blue  edges,  and  edges  meeting 
only  green  faces  are  called  green  edges,  and  vertices  meeting  only  green  faces  are  called  green 
vertices. 

To  require  that  the  blue  faces  have  at  most  4  sides  is  not  essential  but  simplifies  things:  it 
is  baised  on  the  assumption  that  the  blue  faces  are  formed  by  intersecting  sleeve  faces  with  the 
halfspace  above  the  separating  plane,  and  that  5  is  in  general  position  so  all  the  faces  of  H{S) 
are  triangular.  From  our  previous  discussion  we  know  that  by  orienting  the  blue  faces  we  can 
describe  a  circuit  of  directed  edges  which  includes  all  seam  vertices  and  edges  but  can  include 
the  same  (undirected)  edge  twice  (with  both  orientations):  specifically,  when  two  blue  faces  meet 
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Figure  12:  A  seam  polytope:  sea,  mainland,  islands,  peninsulas. 

at  a  seam  edge,  that  edge  will  be  met  twice  in  the  directed  circuit.  Henceforth  the  seam  will 
generally  denote  the  directed  chcmt  so  formed.  It  cannot  be  sissumed  that  every  green  face  meets 
the  top  of  the  seam  polytope;  moreover,  the  green  faces  can  have  any  number  of  bounding  edges, 
and  their  connectivity  pattern  can  be  complex. 

In  this  section  when  we  speak  of  connectivity  we  mean  'face  connectivity.'  In  other  words, 
two  blue  or  green  faces  are  adjacent  in  this  relation  if  and  only  if  they  share  an  edge:  it  is  not 
enough  that  they  have  a  vertex  In  common.  We  speak  of  'green  component'  to  mean  a  connected 
component  of  the  set  of  green  faces  in  this  sense.  The  set  of  blue  faces  is  connected,  since  each 
face  is  connected  to  the  base  (we  visuaUze  it  as  forming  a  connected  blue  'sea');  There  is  one 
green  component  meeting  the  top,  which  will  be  called  the  'mainland.'  however,  there  may  be 
more  green  'islands'  and  'peninsulas'  formed  where  the  seam  returns  upon  itself,  and  the  seam 
may  visit  the  same  edge  twice,  so  both  incident  faces  are  blue  faces.  A  green  component  is  a 
'peninsula'  when  it  shares  a  vertex  with  another  component  at  the  point  where  the  seam  first 
meets  it.  If  at  the  point  where  the  seam  first  meets  it  the  incoming  and  outgoing  seam  edges 
are  coincident  (blue)  edges,  then  we  think  of  the  component  as  an  island.  This  distinction  helps 
us  to  visualize  the  polytope:  it  is  not  important  algorithmicaUy.  See  figure  12. 

Let  us  reconsider  briefly  the  construction  of  the  sleeve  and  the  two  seam  polytopes  that  are 
built  during  this  construction.  Without  loss  of  generality  we  consider  only  the  upper  seam.  The 
base  of  the  seam  polytope  is  the  intersection  of  H(S)  with  the  separating  plane  R.  The  top  of 
the  seam  polytope  may  be  defined  by  a  horizontal  (parallel  to  R)  plane  which  passes  just  above 
the  highest  vertex  on  the  upper  seam  (since  the  vertices  in  S  can  be  assumed  presorted  with 
respect  to  the  i-coordinate,  this  vertex  can  be  obtained  with  little  computational  overhead.) 
The  planes  supporting  the  faces  on  the  sleeve  define  the  blue  faces  of  the  polytope;  the  planes 
supporting  those  faces  of  H{S)  which  meet  the  blue  faces  above  R  define  the  green  faces  of  the 
seam  polytope. 

Note.  It  is  significant  that  the  green  faces  are  faces  of  H{P)  and  the  blue  faces  are  faces 
of  the  sleeve,  but  whereas  a  part  of  every  sleeve  face  becomes  a  blue  face  of  the  polytope,  only 
those  faces  of  H{P)  incident  to  the  seam  from  above  define  green  faces  of  the  seam  polytope. 
Therefore  the  green  faces  will  in  general  be  larger  than  the  corresponding  faces  of  H{P).  The 
advantage  is  that  their  interconnection  structure  is  simpler  than  the  structure  oi  H{P)  itself 


For  the  rest  of  this  section,  if  A'  is  any  polytope  and  A  any  subset  of  K,  suppose  that  we 
have  determined  for  a  given  line  L  that  if  L  penetrates  K  then  it  penetrates  A,  and  if  it  does 
not  penetrate  K  then  the  front  (respectively,  back)  plane  through  L  to  K  must  coincide  with 
the  front  (back)  plane  through  L  to  A.  Then  we  say  that  the  given  line-query  (for  front  or  back 
plaine)  hats  been  localized  to  A.  Since  the  front  plane  to  K  through  L  coincides  with  the  back 
plane  to  K  through  L',  where  L'  is  the  same  line  as  L  with  the  opposite  orientation,  we  can 
always  restrict  our  attention  to  the  front  plane.  (For  our  purposes  A  will  always  be  the  convex 
hull  of  some  of  the  corners  o{  K.) 

A  hierarchy  of  seam  polytopes  can  be  constructed  while  H{S)  is  being  built;  for  any  oriented 
line  L,  a  line-query  for  H{S)  and  L  can  be  solved  by  a  recursive  sequence  of  line-queries  for  L 
and  this  hierarchy  of  seam  polytopes,  Eis  follows. 

Lemma  15  Let  H{S)  be  as  defined  above  and  suppose  that  H  is  the  seam  polytope  for  Hs  upper 
seam.  Given  an  oriented  line  L,  suppose  that  we  wish  to  determine  whether  L  intersects  the 
interior  of  H{S)  and  if  not  to  give  the  front  plane  through  L  tangent  to  H{S).  Then  a  st)lution 
to  the  corresponding  question  for  U  either  solves  the  query  for  H{S)  or  localizes  it  to  H{P)  or 
H[Q). 

Proof.  First  consider  the  cases  where  L  intersects  the  boundary  of  11  at  two  points.  These 
are  either  (a)  through  the  base  and  another  face,  (b)  through  a  blue  face  and  another  face,  or  (c) 
through  a  green  face  and  the  top  or  another  green  face.  Clearly  in  the  first  two  cases  L  intersects 
the  interior  of  H{S)  (the  base  and  the  blue  faces  are  contained  in  H{S)). 

In  case  (c)  let  us  distinguish  three  subcases.  The  first  is  where  L  intersects  some  green  face 
in  some  green  component  A,  not  the  mainland,  and  the  other  point  where  L  meets  the  boundary 
of  n  is  on  another  component  (perhaps  the  mainland).  For  definiteness  assume  yl  is  an  island. 
Imagine  reconstructing  that  part  of  the  boundary  of  H{S)  contained  within  the  same  part  of 
the  seam  as  the  island  A.  This  can  be  done  piecemeal  by  'shaving'  parts  off  the  island  (thus 
introducing  planes  supporting  green  faces  of  H{S)  not  included  in  the  polytope).  When  such  a 
new  face  is  created  it  is  easy  to  see  that  L  continues  to  intersect  the  green  island  (perhaps  in  the 
new  face)  in  exactly  one  point.  Repeating  this  process  we  eventually  obtain  a  polytope  EI'  where 
the  island  A  m  Jl  matches  a  green  island  A'  in  11'  and  the  latter  island  exactly  matches  part 
of  the  (green)  boundary  of  H{S).  Thus  L  intersects  H{S)  in  this  case.  In  the  second  subcase, 
L  intersects  two  green  faces  in  the  same  component  (perhaps  the  mainland).  We  imagine  a 
similar  'shaving'  procedure  by  which  we  eventually  match  the  appropriate  part  of  the  boundary 
of  H{S).  This  time  there  are  two  possibilities:  either  at  the  end  of  the  process  L  is  demonstrated 
to  intersect  H{S)  (in  two  points  in  the  same  component)  or  at  some  step  the  two  intersection 
points  are  separated  from  the  polytope.  At  this  step  the  front  plane  through  L  tangent  to  the 
resulting  polytope  is  tangent  to  the  face  just  introduced:  let  A  be  the  pyramid  just  shaved  off, 
so  L  intersects  it  but  not  its  base  (which  is  also  the  new  face  of  the  resulting  polytope);  the 
polygonal  region  of  intersection  of  A  with  the  given  front  tangent  plane  must  meet  the  shaved 
polytope,  so  it  must  meet  the  new  face.  While  subsequent  'shaving'  operations  may  move  the 
points  of  contact  of  the  corresponding  front  planes,  they  all  maintain  contact  with  the  same 
green  component.  Since  this  represents  a  green  part  of  the  polytope,  it  is  in  H(P)  and  the 
query  has  been  localized  in  both  cases  to  H{P).  In  the  third  subcase,  the  line  L  penetrates  the 
mainland  and  the  top:  this  is  essentially  the  same  as  the  second  subcase. 

Next  consider  the  cases  where  L  does  not  intersect  the  interior  of  11,  and  so  the  front  plane 
T  through  L  tangent  to  11  is  well-defined.  If  T  meets  11  on  the  seam,  then  using  a  point  query 
we  can  determine  whether  it  penetrates  the  interior  o(  H{S)  near  the  point  of  tangency  with  11. 
If  it  does  not  then  it  solves  the  query  for  H{S).  Otherwise,  it  either  (d)  intersects  two  faces  on 
the  sleeve  or  (e)  intersects  two  faces  of  H{P)  above  the  sleeve.  In  case  (d)  L  intersects  H[S)\  in 
Ccise  (e)  either  L  intersects  a  face  of  H{P)  or  the  front  plane  tangent  to  H{P)  solves  the  query 
for  H{S). 

Finally,  the  point  of  tangency  may  be  (f)  on  the  top  or  (g)  on  the  base  in  which  case  the 
query  is  localized  to  H{P)  or  H{Q)  respectively.  Q.E.D. 
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The  following  is  an  immediate  corollary. 

Lemma  16  Suppose  thai  a  hierarchy  of  seam  polyiopes  have  been  constructed  along  wtih  H{S), 
and  that  S  has  size  0(n).  If  auxiliary  search  structures  have  been  provided  which  enable  the  line 
queries  for  the  seam  polyiopes  to  be  solved  m  serial  time  0(log(n))  then  line  queries  for  H{S) 
can  be  solved  in  0(log  (n))  sej^al  time  overall. 

We  shall  next  see  an  appropriate  search  structure  for  an  individual  seam  polytope  which 
achieves  the  time  bound  in  Lemma  15.  If  an  individual  seam  polytope  together  with  a  search 
structure  can  be  constructed  in  CREW  time  0(log"(n))  with  n  processors  —  and  it  can,  3ls  we 
shall  see  —  then  the  overall  time  of  0(log  (n))  for  the  construction  of  H{S)  is  achieved. 

7.2      Nested  structure  of  seam  and  bridges 

In  this  section  we  shall  begin  analyzing  the  structure  of  the  seam  and  the  green  section  of  a 
given  seam  polytope.  We  aim  to  identify  one  or  more  treehke  structures  in  the  set  of  green  edges 
which  will  permit  us  to  decompose  the  polytope  in  a  balanced  way  in  the  spirit  of  Goodrich's 
improvement  to  the  Voronoi  diagram  algorithm  incorporated  into  Section  4.  As  already  noted  the 
polytope  can  be  viewed  as  composed  of  top,  baise,  a  blue  sea,  and  one  or  more  green  components, 
of  which  one,  the  mainland,  surrounds  the  top.  The  secim  can  be  represented  as  a  sequence  of 
vertices,  possibly  with  repetitions.  We  can  cissume  that  this  sequence  is  stored  in  an  array  which 
we  may  caU  the  seam  list.  We  assume  that  the  edge  connecting  the  last  vertex  on  the  list  to  the 
first  is  on  one  of  the  green  faces  which  reach  the  top;  this  implies  that  it  is  on  the  mainland.  The 
notion  of  the  span  of  a  vertex  will  prove  useful:  given  any  seam  vertex  v,  let  i  be  the  first  and 
j  the  last  index  where  it  appears  in  the  seam  Ust.  Then  the  interval  [i . . .  j]  is  called  the  span  of 
V.  Among  the  green  edges  (edges  bounding  two  green  faces)  let  us  call  those  Incident  to  the  top 
or  to  the  seam  connector  edges.  If  we  remove  these  connector  edges,  the  remaining  green  edges 
cannot  contain  any  cycles  (otherwise  there  would  be  green  faces  bounded  away  from  the  seam); 
in  other  words,  they  form  a  forest  of  free  trees. 

Definition  17  A  bridge  is  either  (a)  a  connector  joining  two  seam  vertices  or  (b)  a  connector 
joining  the  top  to  the  seam  or  (c)  one  of  the  free  trees  of  green  edges  and  vertices  together  with  its 
incident  connectors.  Bridges  of  type  (a)  are  called  singular  bridges,  those  of  type  (b)  art  called 
trivial,  and  those  of  type  (c)  art  called  nonsingular. 

See  figure  13.  The  following  lemma  reflects  the  nesting  structure  of  seam  vertices  and  bridges. 
The  first  part  says  that  any  two  seam  vertices  have  disjoint  spans  or  the  span  of  one  includes 
that  of  the  other. 

Lemma  18  Given  distinct  seam  vertices  v  and  v'  with  spans  [i,j]  and  [i' ,  j']  respectively,  if 
i  <  i'  <  j  then  j'  <  j.  Given  any  bridge  B,  either  all  the  seam  vertices  incident  to  B  have 
disjoint  spans  or  this  holds  for  all  but  one  of  these  vertices  where  the  exceptional  vertex  has  a 
span  which  includes  all  the  others. 

Proof.  The  first  part  follows  from  a  simple  topological  argument:  suppose  that  the  seam  visited 
V  at  i  and  j  and  v'  at  t'  and  j'  where  i  <  i'  <  j  <  j' .  Then  one  can  argue  that  the  path  between 
j'  and  j'  includes  within  it  (on  the  'right-hand  side')  at  least  one  blue  face  bounded  away  from 
the  base.  Namely,  the  path  from  i'  to  j  would  have  blue  faces  to  its  immediate  right;  so  would 
the  path  from  j  to  j';  therefore  the  combined  path  would  enclose  some  blue  faces  of  the  polytope: 
since  the  seam  is  bounded  away  from  the  base  these  faces  would  h  '•  ounded  away  from  the  base. 
See  figure  14. 

For  the  second  part,  define  the  depth  of  a  seam  vertex  as  follows:  the  seam  vertices  at  depth 
0  are  those  not  included  in  the  span  of  any  other  vertex;  inductively,  those  at  depth  k  -I-  1  are 
included  in  the  span  of  vertices  of  depth  k  but  none  at  depth  k+  1.  A  seam  vertex  is  at  depth  0  if 
and  only  if  it  is  on  the  mainland.  Vertices  at  greater  depths  are  either  incident  only  to  blue  faces 

29 


Figure  13;  Bridges  on  a  seam  polytope. 


Figure  14:  The  spans  of  seam  vertices  cannot  interlace. 
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Figure  15:  Covering  relation:  1  covers  2,  3  covers  1  and  2,  5  covers  4. 

or  are  on  green  components;  all  vertices  on  the  same  component  are  at  the  same  depth  except 
for  that  vertex  where  the  seam  first  meets  the  component  (and  last  leaves  it).  Since  clearly  all 
seam  vertices  incident  to  the  same  bridge  are  on  the  same  component,  and  all.  seam  vertices  at 
the  same  depth  have  disjoint  spans,  the  rest  of  the  lemma  follows.  Q.E.D. 

Next  we  define  the  inner  vertex  span  of  a  nontrivial  bridge  B.  Let  B  be  a  nontrivial  bridge 
and  V  the  set  of  seam  vertices  on  B:  if  all  vertices  in  V  have  disjoint  spans  then  we  define  the 
leftmost  (rightmost)  seam  vertex  on  B  is  the  vertex  in  V  whose  span  is  leftmost  (respectively, 
rightmost).  On  the  other  hand,  if  some  vertex  in  V  has  a  span  including  the  spans  of  all  others 
in  V  then  it  is  defined  as  the  leftmost  seam  vertex;  among  the  other  vertices  in  V  the  one  whose 
span  is  rightmost  is  defined  as  the  rightmost  seam  vertex  on  B.  The  inner  vertex  span  of  a 
nontrivial  bridge  B  is  the  smallest  interval  [i . . .  j]  where  j  is  the  first  seam  occurrence  of  its 
rightmost  seam  vertex  and  i  the  last  seam  occurrence  before  j  of  its  leftmost  vertex.  Note  that 
the  span  of  the  bridge  B  does  not  necessarily  include  the  spans  of  its  extreme  seam  vertices  but 
it  includes  the  spans  of  all  other  vertices  in  B.  That  part  of  the  seam  represented  by  the  span 
we  call  the  lower  path  for  B.  (Note  that  we  call  this  a  path  for  B\  since  it  follows  the  seam,  it 
is  not  a  path  in  B .)  There  is  a  unique  path  in  B  of  green  edges,  called  the  upper  path  in  B, 
connecting  the  extreme  vertices  such  that  all  vertices  in  B  are  contained  between  the  lower  path 
for  B  and  upper  path  in  B.  (This  will  be  covered  in  more  detail  later.)  Combining  these  two 
paths  (the  upper  being  directed  from  right  to  left)  we  obtain  the  circuit  around  B.  This  circuit 
is  not  necessarily  a  simple  circuit.  The  part  on  it  and  to  its  left  we  consider  the  interior  of  the 
circuit,  which  it  surrounds.  Thus  all  of  S  is  surrounded  by  this  circuit. 

Definition  19  Given  two  distinct  nontnvtal  bridges  B  and  B' ,  B  covers  B'  if  the  circuit  around 
B  surrounds  all  of  B'. 

See  Figure  15.  It  is  clear  from  the  definition  that  this  covering  relation  is  acyclic  and  irreflexive 
and  can  be  considered  to  give  the  set  of  bridges  the  structure  of  '  Tirest  of  rooted  trees.  We 
define  a  principal  bridge  to  be  a  nontrivial  bridge  which  is  on  the  mainland  and  which  is  covered 
by  no  other  bridge.  Thus  a  nonsingular  principal  bridge  is  one  which  is  connected  to  the  top 
of  the  seam  polytope.  There  is  at  most  one  connector  joining  any  bridge  to  the  top,  because 
otherwise  there  would  exist  a  green  face  bounded  away  from  the  seam. 
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Figure  16:  Near-split  of  a  polytope. 

We  give  every  bridge  a  root  as  follows:  in  a  principal  nonsingular  bridge,  the  root  is  the 
unique  vertex  joined  to  the  top.  In  any  other  nonsingular  bridge  the  root  is  on  the  seam:  it 
is  the  leftmost  vertex  on  the  seam  (strictly  speaking  to  avoid  ambiguity  the  root  should  be 
associated  with  the  leftmost  connector  in  clockwise  order  joining  this  vertex  to  the  bridge).  We 
shall  see  later  that  the  covering  relation  among  the  bridges  and  the  parent  relation  within  them 
(implicitly  fixed  by  choosing  roots  for  all  the  bridges  considered  as  trees)  can  be  recovered  from 
information  about  the  order  of  occurrence  of  connectors  along  the  seam. 

7.3     Near  splits,  core  and  lobes  of  a  seam  polytope. 

The  data-structure  we  shall  describe  for  line-queries  will  make  extensive  use  of  a  partition  method 
which  we  call  near-splitUng.  A  natural  divide-and  conquer  method  to  solve  a  line  query  for  a 
line  L  and  a  polytope  A'  would  be  to  divide  A  by  a  plane  into  two  simpler  parts:  by  solving 
the  line-query  relative  to  the  face  common  to  the  two  parts  we  could  localize  to  one  or  other 
part.  Since  such  a  splitting  face  could  have  a  large  number  of  vertices,  to  localize  the  query 
would  require  logarithmic  serial  time.  Therefore  rather  than  sphtting  K  into  two  parts  meeting 
at  a  common  face  we  shall  decompose  K  into  two  convex  polytopes  which  meet  in  a  set  with 
nonempty  interior  and  whose  union  is  K .  This  motivates  the  following  definition. 

Definition  20  Lei  K  be  a  seam  polytope,  and  J  a  polygonal  Jordan  curve  whose  comers  are 
comers  of  K  and  wh  "^e  edges  are  in  ike  boundary  of  K  (so  successive  comers  of  J  are  always 
in  the  same  face).  Lei  A  be  the  convex  hull  of  those  comers  of  K  on  and  outside  J;  let  B  be 
the  convex  hull  of  those  comers  of  K  on  and  inside  J .  Then  A  and  B  form  a  near  spUt  of  the 
polytope  K.  See  figure  16. 

Lemma  21  With  A,  B,  J ,  and  K  as  in  the  above  definition,  let  L  be  an  oriented  line.  Then  if 
suitable  test  polygons  are  provided  for  both  A  and  B  at  each  vertex  of  J ,  a  solution  to  a  line-query 
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Figure  17:  Test  polygons  help  localize  the  query. 

for  L  relative  to  A(^  B  either  solves  the  line-query  for  K  or  can  he  used  with  0{l)  extra  work 
(on  a  serial  processor)  to  localize  the  query  to  A  or  B . 

Proof.  If  one  imagines  A  being  formed  by  successively  discarding  vertices  in  the  boundary 
of  A'  which  are  strictly  inside  J,  it  is  easy  to  conclude  that  J  is  composed  of  edges  and  vertices 
of  A  and  the  boundary  of  A  matches  that  of  K  outside  J  and  that  of  H(J)  inside  J.  Let  ,4i 
(respectively,  .^2)  represent  that  part  of  the  boundary  of  A  on  and  outside  (respectively,  on 
and  inside)  J.  Similarly,  let  5i  (respectively,  B2)  represent  that  part  of  the  boundary  of  B  on 
and  inside  (respectively,  on  and  outside)  J .  Then  Ai  U  Bn  is  the  boundary  of  A'  and  fli  U  A2  is 
simultaneously  the  boundary  of  i/(J)  and  of  ^n5.  Consequently  AuB  —  K  and  AnB  =  H{J). 

Given  a  line-query  for  an  oriented  line  L  against  A,  suppose  that  it  has  been  solved  for  ,4nS. 
If  L  penetrates  this  set  then  clearly  it  penetrates  A':  otherwise,  let  T  be  the  front  tangent  plane 
through  L  to  Ar[B.  Since  this  latter  set  is  H{J),  T  must  be  tangent  to  J:  suppose  that  it  meets 
J  at  some  corner  p.  This  is  on  the  boundary  of  A,  and  since  all  of  J  is  on  the  same  side  of  T,  T 
cannot  cut  both  A  and  B  transversally.  Therefore  to  locaUze  the  query  it  is  enough  to  eliminate 
one  of  these  polytopes  from  consideration:  With  the  aid  of  test  polygons  this  can  be  done  in 
constant  (serial)  time  without  requiring  the  vertex  p  to  have  low  degree.  To  see  this,  let  M  be 
a  plane  which  separates  the  vertex  p  from  the  others  in  both  A  and  B,  so  it  serves  to  define 
test  polygons  for  p  relative  io  Af\  B  (and  to  A  and  B  also).  Consider  the  intersection  of  M 
with  the  boundaries  dA  and  dB  of  A  and  B  respectively.  (See  figure  17).  We  wish  to  determine 
whether  the  line  MnT  meets  dA  -  d{A  n  B)  or  dB  -  d{A  n  B)  transversally  (it  cannot  meet 
both  transversally).  Note  that  the  boundaries  of  A  and  B  meet  at  two  points  c  and  d  of  A/, 
and  the  two  given  curves  lie  on  opposite  sides  of  the  fine  joining  these  two  points.  One  of  these 
curves  can  be  eliminated  in  constant  time  by  considering  the  relative  placements  of  this  hne  and 
the  line  MnT.  Q.E.D. 

N'^w  let  n  be  a  seam  polytope.  We  want  to  supply  it  with  a  structure  suitable  for  solving 
line-<iai;ries:  the  structure  will  be  one  which  uses  'near  splits'  recursively  to  localize  the  queries. 
The  first  step  in  such  a  procedure  will  be  to  decompose  the  polytope  into  a  (rather  simple)  'core' 
with  some  'lobes'  attached. 

Definition  22  A   multiple  vertex  is  a  seam  vertex  visited  more  than  once  by  the  seam.     A 
principal  bridge  is  a  bridge  on  the  mainland  which  is  covered  by  no  other  bridge.     Given  a 
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multiple  vertex  v  there  exist  two  base  vertices  x  and  y  to  its  left  and  right  respectively  and 
adjacent  to  it  (or  at  least  on  blue  faces  incident  to  it)  such  that  the  Jordan  curve  through  vxy 
contains  exactly  thai  pari  of  the  seam  within  the  span  of  v.  We  define  the  lobe  under  v  to  be 
the  subpolyiope  formed  by  all  comers  on  and  mside  this  curve  (thus  containing  the  given  part  of 
the  seam).  Let  B  be  a  singular  principal  bridge  with  two  incident  seam  vertices  u  and  w:  there 
exist  base  vertices  x  and  y  such  that  the  Jordan  curve  uxyw  contains  on  its  inside  exactly  those 
seam  vertices  within  the  span  of  B .  The  convex  hull  of  all  comers  on  and  inside  this  Jordan 
curve  IS  the  lobe  under  B.  Finally  suppose  that  B  is  a  (nonirivial)  nonsingular  principal  bridge, 
which  must  therefore  have  extremal  seam  vertices  u  and  w  and  a  vertex  v  adjacent  to  the  top; 
again  there  exist  base  vertices  x  and  y  on  blue  faces  incident  to  u  and  w  respectively  such  that 
the  Jordan  curve  vuxyw  contains  all  seam  vertices  within  the  span  of  B  on  its  inside;  this  curve 
defines  the  lobe  under  B  in  this  case. 

Given  a  seam  polytope  H,  suppose  that  the  following  operations  are  performed: 

1.  For  each  multiple  vertex  v  on  the  mainland,  a  Jordan  curve  enclosing  the  lobe  under  v  is 
chosen  as  described  above,  and  these  Jordan  curves  are  chosen  not  to  cross  one  another:  let 
n'  be  the  convex  hull  of  all  corners  on  or  outside  aU  these  Jordan  curves.  We  can  consider 
n'  to  have  been  obtained  by  performing  a  lobotomy  at  ea^h  of  these  lobes  —  discarding 
the  vertices  inside  the  given  Jordan  curves.  Each  new  face  of  11'  may  be  considered  a  blue 
face  and  the  resulting  seam  polytope  heis  no  green  components  except  the  mainland. 

2.  For  each  nontrivial  principal  bridge  on  11'  (principal  bridges  on  EI  and  IT'  correspond  one- 
to-one),  choose  a  Jordan  curve  (with  4  or  5  vertices  as  appropriate)  surrounding  the  lobe 
under  that  bridge  and  ensure  that  these  Jordan  curves  do  not  cross  one  another:  let  H" 
be  the  convex  hull  of  all  corners  of  11'  on  or  outside  all  of  these  Jordan  curves.  Thus  the 
lobes  under  each  principal  bridge  have  been  removed. 

3.  For  the  singular  bridges,  let  both  new  faces  created  be  considered  as  blue  faces;  for  the 
nonsingular  bridges,  it  can  be  deduced  from  the  convexity  of  H'  that  uvw  (as  discussed  in 
Definition  22)  forms  a  face,  which  may  be  considered  a  green  face,  and  the  other  two  faces 
may  be  considered  as  blue. 

See  figure  18.  The  polytope  H'  is  called  the  cone  of  11.  This  polytope  has  a  rather  simple 
structure,  since  the  only  green  faces  not  meeting  the  top  are  the  triangular  faces  replacing 
nonsingular  principal  bridges.  All  paths  from  the  top  to  the  base  not  including  seam  edges 
contain  one  blue  edge  and  one  or  two  green  edges.  The  following  lemma  shows  how  to  use 
near-spHts  recursively  to  solve  a  line  query  against  II". 

Lemma  23  Suppose  that  a  seam  polytope  U  is  processed  as  above  yielding  R'  and  then  its  core 
n".  Let  L  be  a  fixed  oriented  line.  Then  in  0(log(n))  serial  time  it  can  be  determined  whether 
L  intersects  the  core,  and  if  not,  a  line-query  for  L  can  be  localized  to  the  correct  lobe  of  the 
polytope  n. 

Proof.  Let  pi  and  po  be  two  vertices  on  the  top,  connected  by  edges  of  the  core  to  vertices  gi 
and  q2,  respectively,  on  the  base,  where  these  connecting  paths  do  not  cross  and  contain  no  seam 
edges  (thus  they  each  have  at  most  three  edges  in  them).  These  can  be  joined  along  the  top 
and  base  to  produce  a  Jordan  curve  which  can  then  be  used  to  form  a  near-spht  of  the  polytope 
into  sub-polytopes  A  and  B.  Although  these  sub-polytopes  are  not  exactly  seam  polytopes, 
they  retain  the  important  feature  that  there  exist  short  paths  (at  most  3  edges)  connecting 
the  base  to  the  top.  Thus  this  procedure  can  be  repeated  in  a  balanced  fashion  to  solve  any 
line-query  against  H"  in  logarithmic  time.  Suppose  that  the  line-query  has  been  solved  in  this 
way,  positively,  in  the  sense  that  a  well-defined  front  plane  T  through  L  tangent  to  11"  at  some 
corner  r  hais  been  found.  Then  x  is  incident  to  at  most  two  lobes  of  11',  and  T  can  meet  the 
interior  of  at  most  one:   hence  the  line-query  can  be  localized  to  the  correct  lobe  of  EI'  using  a 
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Figure  18:  Lobes  in  a  seam  poly  tope. 

point-query.  Next  using  the  structure  defined  later  in  this  section  for  searching  a  lobe,  the  query 
jan  be  solved  against  this  lobe.  Suppose  the  query  yields  a  positive  answer  so  another  plane 
tangent  to  11'  at  some  corner  y  of  11'  is  found.  Then  y  is  incident  to  at  most  one  lobe  of  11  so 
the  same  procedure  can  be  repeated  to  solve  the  query  against  II,  as  required.  Q.E.D. 

7.4      Balanced  decomposition  of  a  lobe. 

We  now  show  how  to  decompose  a  lobe  by  successive  near-splits.  Recall  that  a  lobe  could  have 
been  separated  from  a  seam  polytope  by  a  Jordan  curve  with  3,  4,  or  5  vertices.  The  lobe  has 
essentially  the  structure  itself  of  a  seam  polytope  except  that  it  does  not  reach  the  top  and  it 
has  one,  two,  or  three  'back  faces'  respectively  forming  the  back  of  the  convex  hull  of  the  3  to  5 
vertices  in  question.  These  'back  faces'  do  not  matter  in  line  queries  and  they  need  not  have  any 
particular  color.  Also  the  seam  need  no  longer  be  a  closed  circuit,  but  it  inherits  the  sequencing 
and  nesting  structure  from  the  larger  polytope,  and  the  concept  of  bridge  still  applies.  The 
following  rather  straightforward  lemma  describes  the  basis  of  the  recursive  decomposition. 

Lemma  24  Suppose  that  i  and  j  index  vertices  along  the  seam  of  a  lobe  A  and  P  is  a  path  of 
edges  and  vertices  from  the  vertex  at  j  which  we  shall  denote  v,  say,  to  that  at  i  (denoted  u) 
where  the  comers  of  P  are  (except  at  the  ends)  green  vertices  and  pairs  of  successive  comers 
belong  to  the  same  green  face.  If  u  and  v  are  the  same  vertex  we  assume  thai  P  has  length 
zero.  Then  P  can  be  continued  into  a  Jordan  curve  (with  one  or  two  more  comers  on  it)  in 
the  boundary  of  A  such  that  the  green  veriices  and  edges  surrounded  by  J  are  precisely  those 
surrounded  by  the  seam  path  from  i  to  j  joined  to  P. 

Proof.  Let  e  be  the  seam  edge  (i  —  1,  i)  and  /  the  seam  edge  {j,  j  +  1).  Choose  a  base  point 
X  common  to  the  blue  face  meeting  e;  similarly  choose  a  base  point  y  on  the  same  blue  face  as 
/.  (This  i'  ■>mbiguous  if  and  only  if  e  and  /  are  the  same  edge  with  opposite  orientations,  in 
which  case  we  can  choose  x  and  y  on  the  blue  faces  meeting  this  edge  from  the  left  and  right 
respectively.)  Then  the  path  uxyv  extends  P  to  a  suitable  Jordan  curve:  when  u  and  v  are  the 
same  this  path  is  triangular.  Q.E.D. 

Using  the  hierarchical  structure  among  and  within  bridges  we  can  obtain  a  balanced  decom- 
position of  the  lobe,  but  a  few  details  remain  to  be  considered  first.  The  structure  of  bridges  will 
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Figure  19:  Spanning  face:  bridges,  fictitious  green  edges. 

furnish  paths  P  suitable  for  applying  Lemma  24.  As  already  mentioned  this  structure  wiU  be  ex- 
ploited solely  to  yield  balanced  decompositions  of  a  seam  polytope  or  a  lobe  of  a  seam  polytope. 
In  the  absence  of  multiple  seam  vertices  this  structure  would  furnish  us  with  the  decompositions 
we  need  (after  some  effort):  however,  because  the  seam  can  return  upon  itself  there  could  be 
many  green  components  consisting  of  a  single  green  fau:e  surrounded  by  the  seam  and  incident 
to  no  bridge  whatsoever.  To  answer  this  difficulty  we  shall  introduce  some  'degenerate'  bridges. 
Let  r  be  a  multiple  seam  vertex  and  let  e  and  /  be  the  seam  edges  first  entering  and  last  leaving 
t;.  If  these  edges  are  coincident  or  else  if  the  two  edges  e  and  /  bound  the  same  green  face  on 
their  left,  then  we  regard  the  vertex  v  as  forming  a  degenerate  bridge  whose  span  is  the  span  of 

V. 

A  spanning  face  is  a  green  face  which  the  seam  either  surrounds  completely  or  meets  in  a 
disconnected  set  of  indices. 

In  the  following  lemma,  recall  that  the  bridges  have  been  given  the  structure  of  rooted  trees 
rather  than  just  free  trees  by  making  the  topmost  vertex  of  principal  bridges  and  the  leftmost 
vertex  of  other  bridges  the  root.  (More  precisely,  the  leftmost  connector  on  these  bridges  defines 
the  roots.) 

Lemma  25  Let  F  be  a  spanning  face  in  the  lobe.  Then  (i)  there  is  a  unique  vertex  v  or  bridge 
B  whose  span  includes  all  other  bridges  and  seam  vertices  incident  to  the  face;  (ii)  for  every 
other  bridge  B'  incident  to  F ,  its  first  and  last  seam  vertices  are  incident  to  F,  and  B  is  the 
parent  of  B'  in  the  covering  relation  (i.e.,  B  covers  B'  and  every  other  bridge  covering  B'  also 
covers  B);  (Hi)  if  B  and  B'  are  two  bridges  such  that  B  is  the  parent  of  B'  with  respect  to  the 
covering  relation,  then  either  both  are  degenerate  and  connected  by  a  blue  seam  edge  or  they  have 
a  vertex  in  common  or  there  exists  a  spanning  face  F  such  that  (ii)  holds  for  F,  B,  and  B'. 

Proof,  (i)  The  seam  leaves  the  face  F  as  often  as  it  enters  it:  thus  the  first  and  last  vertices 
where  it  meets  F  either  (a)  coincide  or  (b)  form  the  bounds  on  a  chain  of  one  or  more  green 
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edges  around  F .  In  the  latter  case,  by  definition,  all  edges  in  the  chain  belong  to  the  same  bridge 
B.  This  proves  the  existence  either  of  the  unique  spanning  vertex  in  the  case  (a)  or  the  spanning 
bridge  B  in  case  (b).  (ii)  Now  let  B'  be  any  other  bridge  incident  to  F  in  a  chain  C  of  green 
edges,  say  (or  possibly  at  a  single  seam  vertex).  The  seam  and  the  chain  of  edges  meeting  B  (or 
V  in  the  first  case)  between  them  surround  a  formation  of  one  or  more  green  components.  Each 
green  component  in  this  formation  is  simply  connected;  therefore  C  separates  B'  from  the  rest 
of  this  formation,  and  hence  the  seam  cannot  meet  B'  outside  the  span  of  the  vertices  bounding 
C  (cf.  the  argument  illustrated  in  Figure  14).  This  implies  that  these  are  the  extreme  seam 
vertices  on  B' ,  and  the  leftmost  vertex  in  C  is  the  root  of  the  bridge  B'  under  the  conventions 
previously  established.  Thus  B  covers  B' .  Moreover,  if  B'  covers  some  other  bridge  B"  it  is 
clear  that  B"  is  on  the  other  side  of  the  chain  C  so  B"  cannot  be  incident  to  the  face  F .  (iii)  If 
B  and  B'  are  on  the  same  island  (or  mainland)  or  peninsula,  there  cannot  be  any  chain  of  green 
edges  separating  these  two  bridges  —  otherwise  B  would  not  be  the  psu-ent  of  B':  it  follows 
easily  that  there  exists  a  green  face  to  which  both  bridges  are  incident  and  (iii)  holds  in  this 
case.  If  B  and  B'  are  separated  along  the  seam  by  at  least  one  blue  seam  edge,  then  because 
of  the  existence  of  degenerate  trees  we  are  assured  that  they  are  incident  to  opposite  ends  of  a 
unique  seam  edge.  The  last  case  then  is  where  B'  is  on  a  peninsula  within  the  span  of  B  but  not 
containing  B.  Again  because  of  the  existence  of  degenerate  trees  we  know  that  the  peninsula 
containing  B'  must  have  a  point  in  common  with  the  green  component  containing  S  at  a  point 
X,  say,  and  it  foUows  that  either  x  is  on  5  or  x  is  separated  from  S  by  a  spanning  face  which 
meets  the  seam  within  the  span  of  S.  Q.E.D. 

Now  the  parent  relation  within  the  bridges  themselves  can  be  augmented  to  reflect  the  cov- 
ering relation  among  the  bridges,  as  follows.  Suppose  that  the  vertex  v  or  bridge  B  hcis  been 
defined  as  in  Lemma  25  (ii);  if  it  is  a  bridge  B  rather  than  a  single  vertex  which  is  in  question, 
let  V  be  the  rightmost  seam  vertex  on  the  spanning  face  meeting  B.  Thus  in  either  case  we  have 
distinguished  a  certain  vertex  v.  Then  imagine  the  other  incident  bridges  B'  and  multiple  seam 
vertices  v'  in  anticlockwise  sequence  about  the  face;  connect  the  roots  of  these  other  bridges  to 
t;  by  adding  fictitious  (green)  edges  to  the  lobe,  and  make  v  the  parent  of  all  these  roots.  In  the 
ccise  where  v  has  several  connectors  incident  to  it,  so  it  represents  the  leaf  of  more  than  one  tree, 
ties  may  be  broken  arbitrarily,  so  long  as  only  one  of  the  leaves  corresponding  to  t;  is  chosen  to 
be  parent  and  the  other  leaves  corresponding  to  v  continue  as  leaves  in  the  larger  structure.  See 
figure  19.  Similarly,  the  blue  seam  edges  can  be  incorporated  into  the  structure  by  regarding 
they  themselves  as  forming  links  in  the  tree. 

The  following  lemma  is  then  easily  proved. 

Lemma  26  With  the  addition  of  these  fictitious  edges  and  parent  links,  the  set  of  green  vertices 
and  edges  together  with  the  seam  vertices  suitably  replicated  on  the  lobe  acquires  the  structure  of 
a  single  tree,  whose  root  is  the  topmost  vertex  on  the  lobe. 

Sketch  of  proof.  If  F  is  any  green  face,  then  whether  or  not  it  is  a  spanning  face  we  know  that 
all  bridges  meeting  the  face  are  connected  together  (possibly  through  fictitious  green  edges).  This 
implies  that  all  vertices  on  the  same  green  component  are  connected  together.  The  construction 
ensures  also  that  vertices  in  different  components  are  connected  together  by  chains  of  parent 
links.  Since  the  topmost  vertex  on  the  lobe  is  a  seam  vertex  and  by  definition  it  covers  all  other 
seam  vertices,  it  becomes  the  root  of  the  resulting  binary  structure.  Q.E.D. 

We  are  almost  ready  to  build  a  balanced  decomposition  of  a  lobe.  The  last  requirement  is  to 
embed  the  structure  in  that  of  a  binary  tree.  This  is  quite  Standard:  if  a  node  has  t  >  3  children, 
replace  it  by  ^  -  1  fictitious  nodes,  form  a  complete  binary  tree  with  these  as  internal  nodes 
and  what  were  its  children  as  leaves,  and  let  the  root  of  this  new  tree  correspond  to  the  original 
node.  The  order  of  leaves  in  the  new  tree  matches  the  order  of  the  children  of  the  original  node. 
Let  T  denote  the  resulting  binary  tree  structure.  We  shall  call  the  tree  T  the  binary  structure 
of  the  lobe  for  which  it  was  defined.  Using  the  tree  decomposition  technique  as  described  in 
Section  4,  we  can  find  a  vertex  i;  in  T  such  that  the  subtree  T„  has  between  1/3  and  2/3  as  many 
nodes  as  T.    Let  T^  be  the  tree  T  with  all  proper  descendants  of  T^  removed,  and  recursively 


Figure  20:  Illustrating  cases  (a)  and  (c)  in  Lemma  27. 

decompose  T„  and  T^.  This  can  be  used  to  produce  a  balanced  decomposition  of  the  node  by 
recursive  near-spUts,  basing  the  method  on  the  following 

Lemma  27  Lei  T  be  the  binary  structure  of  a  lobe  A.  IfT  is  decomposed  as  above  into  T„  and 
Tl,  there  exists  a  1-,  2-,  or  3-vertex  near  split  of  A  into  polytopes  A  and  B  such  that  T„  is  the 
binary  structure  of  A  and  T^  thai  of  B. 

Sketch  of  proof.  The  nodes  of  T  correspond  in  a  many-one  fashion  to  the  green  vertices  and 
connectors  on  the  lobe.  If  x  is  a  node  of  T  let  x'  be  the  corresponding  vertex  or  connector  on  the 
lobe.  Given  the  node  v,  let  u  and  w  be  the  first  vertices  encountered  in  following  the  leftmost 
(respectively,  rightmost)  path  down  from  v  such  that  the  corresponding  entity  u'  (respectively, 
w')  represents  a  seam  vertex  (more  ex2u;tly,  a  connector).  It  is  easy  to  see  that  in  the  chain 
P  =  w'v'u'  each  pair  of  successive  vertices  is  in  the  same  green  face  (they  need  not  be  distinct). 
Furthermore  it  is  easy  to  see  that  there  is  an  interval  [i . . .  j]  such  that  the  seam  vertices  indexed 
in  this  interval  consists  precisely  of  those  vertices  x'  such  that  x  is  a  descendant  of  t;  in  T. 
Therefore  by  Lemma  24  there  is  a  Jordan  curve  J  which  surrounds  precisely  those  vertices  (and 
hais  at  most  5  corners).  Suppose  that  B  represents  that  part  within  J  and  A  that  outside  J  in 
a  near-spUt.  Clearly  Tv  is  the  binary  structure  of  B.  The  structure  of  A  is  almost  as  simple: 
it  depends  on  how  many  of  u',  i;',  w'  are  distinct,  (a)  If  they  are  all  distinct  we  have  a  3-vertex 
near  split,  and  u'v'w'  defines  a  triangular  green  face  above  the  seam;  then  T^  almost  reflects 
this  structure:  v  should  be  given  two  children  to  represent  the  new  connectors  in  A.  (b)  If  they 
all  coincide  then  the  Jordan  curve  defines  a  triangular  blue  face  in  A  and  T^  clearly  reflects  the 
structure  of  A.  (c)  If  two  of  them  coincide  (necessarily  u'  =  t;')  then  the  curve  defines  two  new 
blue  faces,  v  represents  a  leaf  vertex  in  the  new  structure  and  again  T^  represents  the  binary 
structure  of  yl.  See  Figure  20.  Q.E.D. 

This  decomposition  proceeds  until  one  arrives  at  a  lobe  with  a  bounded  number  of  green 
faces.  Then  a  line-query  can  be  completed  in  logarithmic  time  because  of  the  following 
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Lemma  28  Lei  A  be  a  lobe  of  a  seam  polytope  with  a  bounded  number  of  green  faces,  and  suppose 
that  a  line-query  has  been  localized  to  this  lobe.  Then  tuith  suitable  preprocessing,  ike  query  can 
be  completed  in  logarithmic  serial  time. 

Proof.  Let  us  consider  the  caise  where  the  polytope  has  no  green  faces:  the  more  general  case  can 
be  covered  with  some  more  attention  to  detail.  In  this  caise  the  polytope  has  a  base,  some  back 
faces,  and  some  blue  faces,  and  all  the  blue  faces  meet  the  base:  the  seam  edges  and  vertices 
now  form  a  treelike  structure  of  blue  edges.  But  this  treeUke  structure  is  subject  directly  to 
Goodrich's  decomposition  technique  and  thus  the  polytope  may  be  further  decomposed  by  the 
methods  already  introduced.  This  ultimately  reduces  the  query  to  a  polytope  with  a  bounded 
number  of  fcices  for  which  the  query  can  be  solved  in  0(1)  serial  time.  Q.E.D. 

7.5      Computing  the  binary  structure  on  a  lobe. 

In  Subsection  7.4  we  saw  that  each  lobe  carried  a  natural  binary  structure  which  could  be 
exploited  to  enable  efficient  solution  of  line-queries  against  the  lobe.  It  is  not  immediately  clear 
how  to  compute  this  binary  structure.  Once  the  interconnection  pattern  among  bridges  has  been 
estabhshed  using  fictitious  green  edges,  to  embed  the  structure  in  a  binary  tree  is  straightforward: 
so  we  concentrate  on  identifying  the  bridges  and  forming  the  fictitious  links. 

To  identify  the  bridges  is  straightforward  once  an  eulerian  cir.cuit  has  been  constructed  to 
traverse  each  bridge  (considered  as  a  free  tree),  and,  moreover,  to  construct  such  a  circuit  can 
easily  be  accomplished  in  0(log(n))  parallel  CREW  time  by  assigning,  say,  one  processor  to  each 
green  vertex  and  one  to  each  green  edge  incident  to  the  vertex;  the  edges  incident  to  the  vertex 
may  be  assumed  sorted  clockwise  to  reflect  the  geometry  of  the  lobe  near  the  vertex.  Once  this 
has  been  done,  all  edges  in  the  same  bridge  can  be  identified  in  logarithmic  parallel  time  using 
a  pointer  doubUng  method.  Thus  we  can  assume  that  each  edge  is  given  a  label  identifying 
the  unique  bridge  containing  it.  Once  this  information  is  available,  it  is  possible  to  identify  the 
correct  root  vertex  on  each  bridge.  Next  assigning  one  processor  to  each  corner  of  each  green  face, 
it  is  possible  to  count  the  total  number  of  seam  vertices  incident  to  each  face:  by  comparing 
this  to  the  range  [i . .  .  j]  representing  where  the  seam  first  and  last  meets  the  face,  one  can 
determine  whether  the  face  is  a  spanning  face.  Given  a  spanning  face  F  it  is  possible  (assigning 
one  processor  to  each  connector  edge  meeting  the  face)  to  indentify  those  bridges  incident  to  the 
face;  and  also  those  'degenerate  bridges'  introduced  above  Lemma  25  can  be  identified  with  one 
processor  per  seam  vertex  meeting  the  face.  All  but  at  most  one  of  these  incident  bridges  will 
have  their  root  vertex  on  the  face,  and  this  information  tells  us  which  vertex  on  the  face  is  to 
be  made  the  parent  of  the  root  vertices  meeting  the  face  using  fictitious  green  edges.  The  order 
of  incidence  of  these  root  vertices  in  anticlockwise  order  around  the  face  can  be  ascertained  in 
logarithmic  parallel  time  using  'parallel  prefix,'  i.e.,  computing  a  cumulative  sum  of  numbers 
around  the  face.  In  this  way  the  interconnections  among  bridges  can  be  established  (to  form 
links  along  blue  seam  edge,  which  the  seam  meets  in  both  directions,  is  a  straightforward  variant 
problem).  It  is  not  difficult  to  splice  the  eulerian  circuits  for  each  bridge  together  with  these 
new  linkages  to  obtain  an  eulerian  circuit  for  the  whole  interconnection  structure.  Finally  the 
correct  parent  relation  for  this  structure  is  easily  computed  using  a  trick  of  Tarjan  and  Vishkin's 
(1985):  break  the  eulerian  circuit  at  the  root  vertex,  then  use  pointer  jumping  to  compute  the 
ranks  of  all  entries  in  the  resulting  linked  list.  Then  each  edge  is  visited  twice  and  it  is  visited 
first  in  the  direction  from  parent  to  child.  We  therefore  conclude  the  following 

Lemma  29  Given  a  lobe  of  a  seam  polytope.  with  the  correct  incidence  relations  established 
among  its  edges,  faces,  and  vertices,  and  with  the  correct  anticlockwise  ordenngs  of  edges  around 
faces  and  vertices,  it  is  possible  to  compute  m  logarithmic  parallel  CREW  time  the  binary  struc- 
ture of  the  lobe. 
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7.6      Building  the  seam  polytope. 

Much  the  same  data-structures  as  have  been  employed  in  the  previous  discussion  will  serve  to 
allow  a  recursive  computation  of  the  seam  polytope  EI.  It  clarifies  things  slightly  if  we  imagine 
the  polytope  to  be  split  in  two  by  a  plane,  so  that  the  seam  is  no  longer  a  closed  circuit:  each 
part  can  be  built  separately  and  they  can  be  combined  at  the  end.  In  the  following  lemma,  the 
irue  seam  of  the  polytope  is  that  part  of  the  seam  which  surrounds  the  mainland. 

Lemma  30  Once  the  seam  of  the  polytope  has  been  found,  the  true  seam  can  be  identified, 
and  the  polytope  can  be  split  into  two  parts  using  a  Jordan  curve  which  divides  the  mainland 
transversally.  This  Jordan  curve  is  defined  by  a  chain  of  four  vertices  connecting  points  on  the 
seam  by  a  path  along  the  base,  together  with  a  plane  whose  intersection  with  the  top  and  the 
mainland  completes  the  circuit. 

Proof.  Once  the  seam  has  been  computed  (and  the  base)  it  is  straightforward  to  construct 
the  true  seam  and  sort  it  in  logarithmic  time  using  a  running-sum  technique.  We  can  therefore 
assume  that  the  polytope  has  only  one  green  component.  Also  a  suitable  plane  supporting  the 
top  can  easUy  be  found,  and  one  green  edge  can  be  identified  on  the  mainland,  namely,  that 
connecting  the  topmost  seam  vertex  p  to  the  top.  Let  C  be  a  plane  through  this  edge  cutting  the 
polytope  transverally  at  the  edge.  Assigning  one  processor  to  each  seam  edge  one  can  identify 
in  logarithmic  parallel  time  the  point  q  closest  to  the  top  where  a  seam  edge  meets  this  plane 
on  the  other  side  of  the  top.  Choose  seam  vertices  x  and  y  sharing  a  blue  face  with  p  and  q 
respectively.  Then  the  plane  C  and  the  chain  pxyq  define  a  suitable  Jordan  curve.  Q.E.D. 

We  can  use  this  Jordan  curve  (formally)  to  split  the  polytope  into  two  overlapping  sub- 
poly  topes  A  and  B.  These  polytopes  have  some  new  'back  fcices'  which  are  not  relevant  to 
the  queries  necessary  for  this  section;  it  is  not  necessEiry  to  know  all  points  along  the  Jordan 
curve  described  in  Lemma  30.  (This  could  be  done  in  0(log  (n))  parallel  time  using  a  suboptimal 
planar  convex  hull  algorithm.  Since  we  wish  to  construct  the  seam  polytope  using  these  structures 
recursively,  this  would  not  be  appropriate.)  Splitting  the  polytope  has  Uttle  importance  except 
that  it  means  that  the  'join'  as  defined  below  is  a  single  zig-zag  Une  like  the  'contour'  occurring 
in  the  construction  of  Voronoi  diagrams  in  Section  4.  First  we  cissume  that  the  green  part  of 
the  polytope  is  connected,  so  the  mainland  is  the  only  green  component.  We  use  divide-and 
conquer:  choose  a  green  seam  edge  e  which  is  midway  in  rank  along  the  seam.  Let  X  be  the  set 
of  green  and  blue  faces  meeting  the  seam  to  the  left  of  this  point  and  Y  the  set  meeting  it  to  the 
right;  recursively  compute  the  polytopes  defined  by  the  top,  base,  back,  and  X  (respectively,  Y), 
and  then  compute  their  intersection.  Let  us  call  the  intermediate  polytopes  in  this  construction 
partial  seam- polytopes.  These  polytopes  are  essentially  seam  polytopes  except  they  may  be 
unbounded  in  one  direction. 

Let  the  join  denote  the  intersection  of  the  'mainlands'  of  two  partial  seam-polytopes.  Clearly 
the  join  meets  the  seam  in  a  single  green  edge  so  it  forms  a  zigzag  chain  of  green  edges  and  plays 
the  role  that  the  contour  played  in  the  context  of  the  2-dimensional  Voronoi  diagram.  Indeed 
the  constructions  are  now  so  similar  that  it  is  enough  to  indicate  some  points  about  the  present 
one. 

Definition  31  A  pure  line-query  is  as  follows:  given  a  semi-infinite  ray  R  known  to  originate 
at  a  point  interior  to  the  core  of  a  partial  seam-polytope  G,  and  not  to  be  contained  within  G,  to 
compute  the  point  at  which  it  leaves  the  polytope. 

Lemma  32  Assuming  thai  the  core  and  lobes  are  available  for  a  partial  seam-polytope  G,  with 
appropriate  balanced  decompositions,  any  pure  line-query  can  be  solved  in  0{\og{n))  serial  time. 
Moreover,  the  exact  coordinates  of  a  solution  vertex  can  be  computed  within  the  stated  time. 

Proof.  Let  R  be  the  ray  given  for  a  pure  line-query.  Testing  R  for  intersection  with  the  top, 
base,  or  any  blue  face,  or  the  back,  can  be  done  in  logarithmic  serial  time  by  binary  search.  If 
so,  the  query  is  solved.  Otherwise,  another  binary  search  will  show  where  R  leaves  the  core;  this 
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solves  the  query  or  localizes  it  to  a  lobe  A.  The  localized  query  can  be  solved  in  logarithmic  time 
35  follows:  suppose  A  =  AuB  is  a  near-split.  If  R  parses  inside  the  Jordan  curve  bounding  AnB 
then  the  query  is  localized  to  the  inner  polytope,  and  moreover,  the  point  where  the  given  ray 
meets  the  back  of  the  inner  polytope  can  be  computed  in  constant  time.  Otherwise  the  query 
is  localized  to  the  outer  polytope;  thus  the  query  is  localized  in  constant  time  by  testing  R  for 
intersection  with  AOB.  Q.E.D. 

Given  X  and  Y  a3  mentioned  above,  let  M  and  TV  be  the  associated  partial  seam-poly  topes 
with  the  appropriate  structures  computed  recursively.  Given  a  vertex  v  of  X,  one  can  find 
whether  d  is  a  vertex  of  the  intersection  as  follows:  let  p  be  some  point  interior  to  the  base  of 
M  n  N:  p  is  easily  computed,  and  serves  as  origin  for  a  family  of  rays  for  pure  line-queries. 

For  any  green  vertex  v  of  A/,  say,  f  is  a  vertex  o{  MnN  if  and  only  if  the  ray  from  p  through 
V  intersects  the  boundary  of  A'^  beyond  v.  Since  the  point  where  the  ray  meets  the  boundary 
can  be  computed  exactly  in  0(log(n))  time  using  the  appropriate  balanced  decompositions,  we 
can  determine  efficiently  whether  or  not  v  belongs  to  the  intersection. 

Thus  the  vertices  of  M  may  be  retained  or  discarded  as  appropriate.  Similarly  for  the  vertices 
of  N .  An  edge  of  M  meets  the  join  if  and  only  if  one  end  v  is  in  the  intersection  and  the  other 
end  w  is  not;  the  point  where  it  meets  the  join  can  then  be  computed  by  solving  a  pure  line-query 
for  the  ray  beginning  at  v  and  passing  through  w.  The  remainder  of  the  'merge'  construction  is 
very  similar  to  that  already  described  for  the  Voronoi  diagram  and  we  omit  the  details. 

This  covers  how  to  construct  the  seam  polytope  under  the  assumption  that  there  are  no 
green  components  except  the  mainland.  Notice  that  such  components  are  eiscertainable  from 
the  structure  of  the  seam  without  constructing  the  seam  polytope,  as  is  the  true  seam,  and  it 
follows  that  a  superset  of  the  seam  polytope  (that  obtained  by  deleting  those  green  faces  not  on 
the  mainland)  can  be  constructed  in  0(log  (n))  parallel  time.  This  polytope  possesses  all  the 
blue  faces  which  the  final  seam  polytope  will  possess. 

So  far  we  have  seen  how  to  compute  the  seam  polytope  corresponding  to  the  true  seam.  It 
remains  to  account  for  the  presence  of  other  green  components.  To  do  this  we  would  proceed 
as  follows.  For  every  green  component  separate  from  the  mainland,  one  can  find  its  topmost 
vertex  v  and  two  base  vertices  u  and  w  which  are  adjacent  to  t;  (i.e.,  which  share  a  face  with 
v)  and  such  that  the  Jordan  curve  vuw  surrounds  all  those  green  components  within  the  span 
of  V.  Thus  using  near-splits  judiciously  one  can  produce  a  formation  of  lobes,  each  of  which 
possesses  just  one  green  component.  Within  each  lobe  the  green  component  can  be  built  using 
the  methods  outlined  above  for  constructing  the  mainland.  (Again  it  is  convenient  to  split  each 
structure  into  two  sections  using  the  method  described  in  Lemma  30.  It  is  not  difficult  to  achieve 
this  using  the  same  methods  as  given  in  the  Lemma.)  Since  this  fragmentation  of  parts  of  the 
seam  polytope  increases  the  total  number  of  vertices  in  all  the  lobes  by  at  most  a  linear  factor, 
the  construction  can  be  accomplished  in  the  same  parallel  time  with  n  processors  overall. 

Thus  we  have  an  0(log"(n))  CREW  construction  of  the  seam  polytope  which  leads  to  the 
C>(log  (n))  construction  of  the  convex  hull,  as  asserted. 

7.7     A  simpler  approach. 

This  concludes  our  discussion  of  the  convex  hull  problem.  The  construction  given  here,  while 
elaborate,  matches  the  time-bounds  of  Anita  Chow's,  and  does  not  require  an  optimal  parallel 
sorting  algorithm.  (Recall  that  the  convex  hull  algorithm  in  [Chow  (1980)]  takes  0(log^(n)  loglog(n)) 
parallel  time,  but  the  factor  of  loglog(n)  arises  from  sorting  and  it  can  be  removed  by  invok- 
ing optimal  parallel  sorting  algorithms.)  We  used  information  about  the  seams  involved  in  the 
recursive  construction  to  provide  a  search-str  .cure  to  answer  hne-queries,  and  we  introduced 
'near-splitting'  as  a  means  to  localize  queries  on  our  seam  polytopes.  The  method  of  near- 
splitting  applies  to  any  polytope  whose  corners  can  be  evenly  partitioned  by  a  short  polygonal 
Jordan  curve  around  its  boundary. 

The  algorithm  presented  in  the  FOCS  paper,  worse  by  a  factor  of  log(n),  depended  on  adapt- 
ing a  variant  of  Kirkpatrick's  planar  decomposition  method  to  a  restricted  kind  of  planar  graph. 
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Richard  Cole  (private  communication)  has  suggested  a  simple  method  to  apply  Kirkpatrick's  de- 
composition to  any  planar  graph  in  parallel.  Let  us  say  only  that  the  decomposition  wais  bcised 
on  the  observation  that  every  planar  graph  has  a  substantial  fraction  of  low-degree  vertices  (ver- 
tices of  degree  k  <  11,  say:  at  least  half  the  vertices  have  low  degree  in  this  sense).  To  develop 
a  suitable  search-structure,  compute  a  maximal  independent  set  of  low-degree  vertices,  discard 
them,  then  triangulate  and  recursively  reprocess  the  resulting  planar  graph.  For  our  application 
the  'planar  graph'  in  question  would  be  dual  to  a  given  convex  polyhedron  K  (see  [Dobkin  and 
Kirkpatrick  (1983)]).  Richard  Cole's  suggestion  is  as  follows:  to  compute  a  maximal  independent 
set  of  low-degree  vertices,  let  G'  be  the  graph  formed  by  the  edges  incident  to  the  low-degree 
vertices.  Construct  a  spanning  forest  of  G',  and  let  S  be  the  subset  of  odd-depth  vertices  in 
this  spajining  forest  (relative  to  an  arbitrary  sissignment  of  roots  for  the  trees  in  this  forest). 
Since  the  trees  have  bounded  degree,  the  size  of  5  is  a  fixed  fraction  of  the  size  of  G' .  Let  G" 
be  the  subgraph  of  G'  (or  equivalently  G)  spanned  by  S;  then  every  vertex  in  G"  has  degree 
at  most  k  —  2  (i.e.,  9).  Repeating  this  process  at  most  five  more  times  we  obtain  a  maximal 
independent  set.  This  method  can  be  implemented  to  solve  the  convex  hull  problem  with  n 
proceessors  in  parallel  time  log"(n)a(n)  where  a{n)  is  the  cost  of  constructing  a  spanning  forest: 
for  instance,  Shiloach  and  Vishkin's  (1982)  CRCW  algorithm  may  be  used  in  which  case  a(n) 
=  log(n),  thereby  yielding  a  much  simpler  NC^{n)  algorithm  for  the  3-D  convex  hull  problem, 
though  on  a  somewhat  more  powerful  model  of  parallehsm  (concurrent  write;  the  policy  by  which 
concurrent  writes  are  resolved  is  immaterial  to  Shiloach  and  Vishkin's  algorithm). 

Recently  Dadoun  and  Kirkpatrick  (1987)  showed  that  a  Kirkpatrick  decomposition  could  be 
built  in  0(log(n)  log'(n))  CREW  time,  thereby  constructing  the  convex  hull  with  an  elegant 
0(log^(n)  log"  (n))  CREW  algorithm. 
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8.  Final  Remarks. 


Although  parallel  algorithms  in  areas  such  as  graph-theoretic  algorithms  and  numerical  problems, 
have  been  well-studied  for  several  years  already,  a  similar  development  is  only  beginning  in  computational 
geometry.  Consequently,  even  the  most  basic  problems  such  as  triangulation  problems  and  those  in 
geometric  optimization  have  not  been  put  in  the  class  NC  until  this  paper.  Standard  techniques  in  the  sub- 
ject such  as  contour-tracing,  plane-sweeping  and  gift-wrapping  initially  seem  inherently  sequential;  one  of 
our  main  contribution  is  to  show  that  //C-analogues  of  these  techniques  actually  exist  As  hinted  before, 
efficient  solutions  to  the  basic  problems  imply  correspondingly  efficient  solutions  to  a  multitude  of  other 
problems,  and  we  have  only  shown  some  of  the  reductions. 

It  is  interesting  to  note  that  in  the  sequential  setting,  the  problems  of  planar  convex  hulls,  planar 
Voronoi  diagrams,  and  3-dimensional  convex  hulls  all  have  the  same  &{n  log/i)  time  complexity.  In  the 
parallel  case,  these  problems  have  been  shown  to  be  in 

NCUn).  NC*2(n),  A^C^(n) 
respectively.    Our  planar  convex  hull  algorithm,  culminating  several  previous  works,  turns  out  to  be 

optimal  in  a  very  strong  sense  and  it  remains  open  whether  the  parallel  solutions  for  the  other  problems  can 
be  further  improved. 

Although  all  common  problems  in  computational  geometry  are  in  NC  [Kozen  and  Yap  (1985)],  we 
think  that  a  main  goal  of  parallel  studies  in  computational  geometry  is  to  bring  problems  into  the  practical 
subclasses  of  NC.  In  this  paper  we  focused  of  NC\in)  for  small  k.  It  is  our  belief  that  such  algorithms  (as 
opposed  to  NC\{n^),  say)  will  have  practical  implications.  In  particular,  it  is  possible  that  they  could  be 
employed  in  VLSI  and  in  graphics  machines  or  robots. 
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